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INTRODUCTION 


This  document  is  the  first  draft  of  the  final  report  for  the  program  New  Methods  in  Robust  Control.  The 
emphasis  of  this  program  was  to  develop  mathematical  theory  to  help  control  system  designers  faced  with  chal¬ 
lenging  control  problems  associated  with  advanced  aerospace  vehicles.  Relevant  applications  include  flight  con¬ 
trol  systems  for  new  Air  Force  fighter/bomber  aircraft,  the  F- 18  HARV  research  vehicle,  the  NASP  vehicle,  the 
next  generation  launch  system  (ALS  or  NLS),  and  the  Space  Station.  A  common  set  of  feature:  characterizing 
these  problems  are: 

1)  Operation  of  vehicles  in  extreme  ranges  of  flight  environment 

2)  Significant  uncertainties  in  dynamic  models  (e.g.  aerodynamics) 

3)  Wide  range  of  parameter  variation  during  flight  (e.g.  mass  properties) 

4)  Performance  driven  system  designs  (small  safety  margins) 

5)  Stability  of  system  is  critical  to  avoid  loss  of  vehicle  and  human  life. 

Robustness  of  a  control  system  is  defined  to  be  its  capability  to  provide  adequate  performance  in  the  presence  of 
uncertainty.  The  importance  of  robusmess  in  aerospace  systems  is  well  appreciated  by  the  people  who  have  to 
fly  them.  Bill  Dana,  NASA  Dryden  Chief  Test  Pilot,  recently  gave  a  presentation  to  the  NASP  community  con¬ 
cerning  lessons  learned  in  the  X-15  Flight  Test  Program  as  they  would  apply  to  NASP.  His  two  key  messages 
were; 

1)  Make  the  vehicle  and  its  controls  very  robust  with  significant  performance  margins 

2)  Make  very  small  incremental  steps  in  the  flight  envelt^  ("0.5  Mach)  during 
the  flight  test  program. 

The  theory  developed  on  our  program  is  motivated  mainly  by  the  need  for  robust  controls  in  teal  vehicles. 

Consistent  with  the  goal  of  this  program,  we  present  three  topics  out  of  those  worked  on  during  the  three  years 
of  this  contract  The  three  topics  are: 

1)  IC-optimal  Control  Theory 

2)  Structured  Singular  Values 

3)  Dynamic  Inversion  Control 

The  first  of  these  topics  is  presented  in  a  paper  written  by  John  Doyle  and  Keith  Glover.  The  H..  theory 
presented  there  is  a  culmination  of  the  research  in  that  area  that  has  been  going  on  for  the  last  decade.  The 
theory  in  its  "final"  form  is  elegant  from  a  mathematical  viewpoint,  but  the  practical  value  of  that  research  lies 
beyond  the  solution  of  the  H_-optimization  problem  itself  (which  has  no  guaranteed  robustness  properties).  H„- 
optimization  is  one  of  the  primary  ingredients  of  the  structured  singular  value  technique.  The  primary  practical 
value  of  the  H,.  theory  is  that  we  can  now  perform  structured  singular  value  synthesis  more  readily. 

The  second  topic  is  presented  in  a  document  written  by  Mike  Elgersma  on  mapping  the  mass-properties  varia¬ 
tions  on  the  Space  Station  into  a  perturbation  structure  for  a  structured  singular  value  design.  Elgersma’s  effort 
to  construct  the  perturbation  structure  was  paid  for  by  Space  Station  contract  funds  -  our  contract  only  contri¬ 
buted  the  extra  funCs  required  for  him  to  document  the  construction  in  a  presentable  form.  'Fhere  are  now 
several  papers  in  the  (pen  literature  on  structured  singular  value  theory  but  there  are  far  too  few  showing  how 
the  perturbation  structure  is  made.  We  included  Elgersma’s  example  here  as  one  application  wlicrc  the  full 
power  of  the  mathematical  theory  can  be  applied  to  a  real-w(M-ld  (real-space?)  application. 

The  third  topic  is  presented  in  ?  Hnrnmcnt  .«epared  by  P'ais"  Merton  xed  Dale  tinns  on  die  subject  of  dynamic 
inversion  for  aircraft  pitch-axis  control.  An  attract  of  this  paper  was  prepared  for  the  Washington  University/ 
AFOSR  woilcshop  in  St.  Louis,  August  15  and  16,  1991.  The  main  new  feauire  of  this  work  is  the  global  stabil¬ 
ity  results  for  complementary  dynamics.  This  result  is  perhaps  the  first  non-lcxal  control  stability  result  that 
applies  directly  to  the  nonlinear  models  used  in  industry  for  mcxlem  aircraft  design. 
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Abstract 

Simple  state-space  formulae  are  derived  for  all  controllers  solving  a  standard  Hoo 
problem:  for  a  given  number  7  >  0,  find  all  controllers  such  that  the  Hoo  norm  of  the 
closed-loop  transfer  function  is  <  7.  Under  these  conditions,  a  parametrization  of  all 
controllers  solving  the  problem  is  given  as  a  linear  fractional  transformation  (LPT)  on  a 
contractive,  stable  free  parameter.  The  state  dimension  of  the  coefficient  matrix  for  the 
LPT  equals  that  of  the  plant,  and  has  a  separation  structure  reminiscent  of  classical  LQG 
(i.e.,  Hi)  theory.  Indeed,  the  whole  development  is  very  reminiscent  of  earlier  Hi  results, 
especiedly  those  of  Willems  (1971).  This  paper  directly  generalizes  the  results  in  Doyle, 
Glover,  Khargonekar,  and  Prancis,  1989,  and  Glover  and  Doyle,  1988.  Some  aspects  of 
the  optimal  case  (<  7)  are  considered. 


1  Introduction 

1.1  Overview 

The  'Hen  norm  defined  in  the  frequency-domain  for  a  stable  transfer  matrix  G(s)  is 
||G||oo  :=  sup  <7[G(ja;)]  (  ct  :=  maximum  singular  value  ) 

U/ 

The  problem  of  analysis  and  synthesis  of  control  systems  using  this  norm  arises  in  a  number 
of  ways.  We  assume  the  reader  either  is  familiar  with  the  engineering  motivation  for  these 
problems,  or  is  interested  in  the  results  of  this  paper  for  some  other  reason.  This  paper 
considers  particulsu-  Tioo  optimal  control  problems  that  are  direct  generalizations  of  those 
considered  in  Doyle,  Glover,  Khargonekar,  and  Francis  (1989),  and  Glover  and  Doyle  (1988), 
hereafter  referred  to  as  DGKF  and  GD,  respectively. 

The  bacic  block  diagram  used  in  this  paper  is 

'This  article  appeared  ip  Dcczdc'  cf  Irfzthcruiiitical  oysiems  ,  A  Ouli€i^tioTi  of  SuTve^^s  uL  iitc 

Occasion  oj  Ute  iOth  birthday  of  Jan  C.  Willems,  H.  Nijmeijer  and  J.M.  Schumacher  (Eds.),  Springer- Verlag 
Lecture  Notes  in  Control  and  Information  Sciences  vol.  135,  1989. 
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where  G  is  the  generalized  plant  and  K  is  the  controller.  Only  finite  dimensional  linear  time- 
invariant  (LTI)  systems  and  controllers  will  be  considered  in  this  paper.  The  generalized 
plant  G  contains  what  is  usually  called  the  plant  in  a  control  problem  plus  all  weighting 
functions.  The  signal  w  contains  all  external  inputs,  including  disturbances,  sensor  noise,  and 
commands,  the  output  z  is  an  error  signal,  y  is  the  measured  variables,  and  u  is  the  control 
input.  The  diagram  is  also  referred  to  as  a  linear  fractional  transformation  (LFT)  on  K  and 
G  is  called  the  coefficient  matrix  for  the  LFT.  The  resulting  closed  loop  transfer  function 
from  w  to  z  is  denoted  by  T,u,  =  ^t{G,  K)  . 

The  main  'Hoo  output  feedback  results  of  this  paper  as  described  in  the  abstract  are 
presented  in  Section  4.  The  proofs  of  these  results  exploit  the  “separation”  structure  of  the 
controller.  If  full  information  (x  and  w)  is  available,  then  the  central  controller  is  simply  a 
gain  matrix  Foo,  obtained  through  finding  a  certain  stable  invariant  subspace  of  a  Hamiltonian 
matrix.  Also,  the  optimal  output  estimator  is  an  observer  whose  gain  is  obtained  in  a  similar 
way  from  a  dual  Hamiltonian  matrix.  These  special  cases  are  described  in  Section  3.  In  the 
general  output  feedback  case  the  controller  can  be  interpreted  as  an  optimal  estimator  for 
FooX.  Furthermore,  the  two  Hamiltonians  involved  in  this  solution  can  be  associated  with  full 
information  and  output  estimation  problems. 

The  proofs  of  these  results  are  constructed  out  of  a  series  of  lemmas,  several  of  which 
have  some  independent  interest,  particularly  those  involving  state-space  characterizations  of 
mixed  Hankel-Toeplitz  operators.  A  possible  contribution  of  this  paper,  beyond  the  new 
formulae  and  theorems,  may  be  some  of  this  technical  machinery,  most  of  which  is  developed 
in  Section  2.  The  result  is  that  the  proofs  of  both  the  theorems  and  the  lemmas  leading  to 
them  «ire  quite  short.  Furthermore,  the  development  is  reasonably  self-contained,  and  the 
primary  background  required  is  a  knowledge  of  elementary  aspects  of  state-space  theory,  £2 
spaces,  and  operators  on  £2,  including  projections  and  adjoints.  More  specialized  knowledge 
about  the  connections  between  Riccati  equations,  spectral  factorization,  and  Hamiltonian 
matrices  would  also  be  useful. 

As  mentioned,  this  paper  is  a  direct  generalization  of  DGKF,  and  contains  a  substantial 
repetition  of  material.  Roughly  speaking,  we  prove  those  results  in  GD  which  were  stated 
without  proof,  using  DGKF  machinery,  which  considered  a  less  general  problem.  An  alterna¬ 
tive  approach  in  relaxing  some  of  the  assumptions  in  DGKF  is  to  use  loop-shifting  techniques 
as  in  Zhou  and  Khargonekar  (1988),  GD,  and  more  completely  in  Safonov  et  al.  (1989).  We 
also  organize  this  paper  much  differently  than  DGKF.  The  results  are  presented  in  a  conven¬ 
tional  bottom-up  linear  order,  with  lemmas  and  theorems  followed  by  their  proofs,  which  in 
turn  only  use  lemmcis  and  theorems  already  proven.  Readers  interested  in  pursuing  all  the 
details  of  the  proofs  may  find  it  more  convenient  than  DGKF.  This  paper  lacks  the  tutorial 
flavor  of  DGKF  and  the  explicit  connections  with  the  more  familiar  H2  problem,  although 
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the  'H2  theory  will  be  found  lurking  at  every  corner. 

We  also  consider  some  aspects  of  generalizations  to  the  <  case,  primarily  to  indicate 
the  problems  encountered  in  the  optimal  ca^e.  A  detailed  derivation  of  the  necessity  the 
generalized  conditions  for  the  Full  Information  problem  is  given.  In  keeping  with  the  style 
of  GD  and  DGKF,  we  don’t  present  a  complete  treatment  of  the  <  case,  but  leave  it  for 
yet  another  day.  Complete  derivations  of  the  optimal  output  feedback  case  can  be  found  in 
Glover  et  al.  (1989)  using  different  techniques. 

1.2  Historical  perspective 

This  section  is  not  intended  as  a  review  of  the  literature  in  "Hoo  theory,  nor  even  an  attempt 
to  outline  the  work  that  most  closely  touches  on  this  paper.  For  a  bibliography  and  review 
of  the  early  'Htx>  literature,  the  interested  reader  might  see  [Francis,  1987]  and  [Francis  and 
Doyle,  1987],  and  an  historical  account  of  the  results  leading  up  to  those  in  this  paper  may  be 
found  in  DGKF.  Instead,  we  will  offer  a  slightly  revisionist  history,  which  lacks  some  factual 
accuracy,  but  has  the  advantage  of  more  clearly  emphasizing  state-space  methods  and,  more 
specifically,  Willems’  central  role  in  'Hoo  theory.  This  mildly  fictionalized  reconstruction  tells 
things  as  they  could  have  been,  if  only  we’d  been  more  clever,  and  thus  contains  a  certain  truth 
as  valuable  as  that  of  a  more  factually  accurate  accounting.  Besides,  “historical  perspectives” 
are  often  revisionist  anyway,  we’re  just  admitting  to  it. 

Zames’  (1981)  original  formulation  of  Hoo  optimal  control  theory  was  in  an  input-output 
setting.  Most  solution  techniques  available  at  that  time  involved  analytic  functions  (Nevanlinna- 
Pick  interpolation)  or  operator-theoretic  methods  [Sarason,  1967;  Adamjan  et  ai,  1978;  BaD 
and  Helton,  1983].  Indeed,  Hoo  theory  seemed  to  many  to  signal  the  beginning  of  the  end 
for  the  state-space  methods  which  had  dominated  control  for  the  previous  20  years.  Unfor¬ 
tunately,  the  standard  frequency-domain  approaches  to  Hoo  started  running  into  significant 
obstacles  in  dealing  with  multi-input-output  (MIMO)  systems,  both  mathematically  and  com¬ 
putationally,  much  as  the  H2  theory  of  the  1950’s  had. 

Not  surprisingly,  the  first  solution  to  a  general  rational  MIMO  Hoo  optimal  control  prob¬ 
lem,  presented  in  [Doyle,  1984],  relied  heavily  on  state-space  methods,  although  more  as  a 
computational  tool  than  in  any  essential  way.  The  steps  in  this  solution  were  as  follows; 
parametrize  all  internally-stabilizing  controllers  via  [Youla  et  ai,  1976];  obtain  realizations  of 
the  closed-loop  transfer  matrix;  convert  the  resulting  model-matching  problem  into  the  equiv¬ 
alent  2  X  2-block  general  distance  or  best  approximation  problem  involving  mixed  Hankel- 
Toeplitz  operators;  reduce  to  the  Nehari  problem  (Hankel  only);  solve  the  Nehari  problem 
by  the  procedure  of  Glover  (1984).  Both  [Francis,  1987]  and  [Francis  and  Doyle,  1987]  give 
expositions  of  this  approach,  which  will  be  referred  to  as  the  “1984”  approach. 

In  a  mathematical  sense,  the  1984  procedure  “solved”  the  Hoo  optimal  control  problem. 
Unfortunately,  it  involved  a  peculiar  patchwork  of  techniques  and  the  associated  complexity 
of  computation  was  substantial,  involving  several  Riccati  equations  of  increasing  dimension, 
and  formulae  for  the  resulting  controllers  tended  to  be  very  complicated  and  have  high  state 
dimension.  Nevertheless,  much  of  the  subsequent  work  in  Hoo  control  theory  focused  on  the 
2  X  2-block  problems,  either  in  the  model-matching  or  general  distance  forms.  This  continued 
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to  provide  a  context  for  a  stimulating  interchange  with  operator  theory,  the  benefits  of  which 
will  hopefully  continue  to  accrue.  But  from  a  control  perspective,  the  Hoo  theory  seemed 
once  again  to  be  headed  into  a  cul-de-sac,  but  now  with  a  in  the  corner.  The  solution  has 
turned  out  to  involve  an  even  more  radical  emphasis  on  state-space  theory. 

In  addition  to  providing  controller  formulae  that  are  simple  and  expressed  in  terms  of 
plant  data,  the  methods  in  DGKF  and  this  paper  are  a  fundamental  departure  from  the 
earlier  work  described  above.  In  particular,  the  Youla  parametrization  and  the  resulting 
2  X  2- block  model- matching  problem  of  the  1984  solution  are  avoided  entirely,  replaced  by 
a  more  purely  state-space  approach  involving  observer- based  compensators,  a  pair  of  2  x  1 
block  problems,  and  a  separation  argument.  The  operator  theory  still  plays  a  central  role  (as 
does  RedheflFer’s  work  [Redheffer,  1960]  on  linear  fractional  transformations),  but  its  use  is 
more  straightforward.  The  key  to  this  was  a  return  to  simple  and  familiar  state-space  tools, 
in  the  style  of  Willems  (1971),  such  as  completing  the  square,  and  the  connection  between 
frequency  domain  inequalities  (  e.g.  11Cj||oo  <  1),  Riccati  equations,  and  spectral  factorization. 
In  essence,  one  only  needed  to  think  about  how  Willems  would  do  it,  and  the  rest  is  simply 
technical  detail. 

The  state-space  theory  of  Hoo  can  be  carried  much  further,  by  generalizing  time-invariant 
to  time-varying,  infinite  horizon  to  finite  horizon,  and  finite  dimensional  to  infinite  dimen¬ 
sional.  A  flourish  of  activity  has  begun  on  these  problems  and  the  already  numerous  results 
indicate,  not  surprisingly,  that  many  of  the  results  of  this  paper  generalize  mutatis  mutandis, 
to  these  cases.  In  fact,  a  cynic  might  express  a  sense  of  dijd.  vu,  that  despite  all  the  rhetoric, 
Hoo  theory  has  come  to  look  much  like  LQG,  circa  1970  (or  even  more  specifically,  LQ  differ¬ 
ential  games).  A  more  charitable  view  might  be  that  current  Hoo  theory,  rather  than  ending 
the  reign  of  state-space,  reaffirms  the  power  of  its  computational  machinery  and  the  wisdom 
of  its  visionaries,  exemplified  by  Jan  Willems. 

1.3  Notation 

The  notation  is  fairly  standard.  The  Hardy  spaces  H2  and  consist  of  square-integrable 
functions  on  the  imaginary  axis  with  analytic  continuation  into,  respectively,  the  right  and  left 
half-plane.  The  Hardy  space  Hoo  consists  of  bounded  functions  with  analytic  continuation 
into  the  right  half-plane.  The  Lebesgue  spaces  £2  =  £2(-oo»oo),  £2-1-  =  £2(0,00)  ,  and 
£2-  =  £2(— oo,0]  consist,  respectively  of  square-integrable  functions  on  (-00,00),  [0,oo), 
and  (—00,0],  and  £„  consists  of  bounded  functions  on  (-00,00).  As  interpreted  in  this 
paper.  Coo  will  consist  of  functions  of  frequency,  C2+  and  C2-  functions  of  time,  and  £2  will 
be  used  for  both. 

We  will  make  liberal  use  of  the  Hilbert  space  isomorphism,  via  the  Laplace  transform  and 
the  Paley- Wiener  theorem,  of  £2  =  £2-1-  ©  £2-  time-domain  with  C2  =  H2  @  H2  in 

the  frequency-domain  and  of  £2+  with  H2  and  £2-  with  Tfj  •  ^  we  will  normally  not 
make  any  distinction  between  a  time-domain  signal  and  its  transform.  Thus  we  may  write 
w  6  £2+  and  then  treat  in  as  if  in  G  7^2 •  This  style  streamlines  the  development,  as  well  as 
the  notation,  but  when  any  possibility  of  confusion  could  arise,  we  will  make  it  clear  whether 
we  are  working  in  the  time-  or  frequency-  domain. 
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For  a  matrix  M  €  M'  denotes  its  conjugate  transpose,  a{M)  =  denotes  its 

maximum  singular  value,  p{M)  denotes  its  spectral  radius  (if  p  =  r),  and  Afl  denotes  the 
Moore-Penrose  pseudo-inverse  of  M.  Im  denotes  image,  ker  denotes  kernel,  and  G'"{s)  := 
G(— J)'.  For  operators,  F*  denotes  the  adjoint  of  F.  The  prefix  B  denotes  the  open  unit  ball 
and  the  prefix  Ttc  denotes  complex-rational. 

The  orthogonal  projections  P+  and  P_  map  £2  to,  respectively,  and  (or  £2+  and 
£2_)-  For  G  G  £»)  the  Laurent  or  multiplication  operator  Me  :  £2  — »  £2  for  frequency- 
domain  m  G  £2  is  defined  by  Mow  =  Gw.  The  norms  on  Coo  and  £2  in  the  frequency- 
domain  were  defined  in  Section  1.1.  Note  that  both  norms  apply  to  matrix  or  vector-valued 
functions.  The  unsubscripted  norm  ||  •  ||  will  denote  the  standard  Euclidean  norm  on  vectors. 
We  will  omit  zJl  vector  and  matrix  dimensions  throughout,  and  assume  that  aU  quantities 
have  compatible  dimensions. 


1.4  Problem  statement 

Consider  the  system  described  by  the  block  diagram 


Both  G  and  K  are  complex-rational  and  proper,  K  is  constrained  to  provide  internal  stability. 
We  will  denote  the  transfer  functions  from  u;  to  z  as  in  general  and  for  a  feedback  connec¬ 
tion  (LFT)  as  above  we  also  write  T,t„  =  Pt{G,  K)  .  This  section  discusses  the  assumptions 
on  G  that  will  be  used.  In  our  application  we  shall  have  state  models  of  G  and  K.  Then 
interned  stability  will  mean  that  the  states  of  G  and  K  go  to  zero  from  all  initial  values  when 
u;  =  0. 

Since  we  will  restrict  our  attention  exclusively  to  proper,  complex-rational  controllers 
which  are  stabilizable  and  detectable,  these  properties  will  be  assumed  throughout.  Thus  the 
term  controller  will  be  taken  to  mean  a  controller  which  satisfies  these  properties.  Controllers 
that  have  the  additional  property  of  being  internally-stabilizing  will  be  said  to  be  admissible. 
Although  we  are  taking  everything  to  be  complex,  in  the  special  case  where  the  original  data 
is  real  (e.g.  G  is  real-rational)  then  all  the  of  the  results  (such  as  K)  will  also  be  real. 

The  problem  to  be  considered  is  to  find  all  admissible  K{s)  such  that  ||Ti„,||oo  <  7  (<  7). 
The  realization  of  the  transfer  matrix  G  is  taken  to  be  of  the  form 
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G(s)  = 


■  A 

Bi 

B2 

Cl 

Du 

Bi2 

C2 

B21 

0 

'  A 

B  ' 

C 

B 

compatible  with  the  dimensions  z(t)  6  ,  y(t)  6  w(t)  E  ,  u(t)  E  and  the  state 

x(t)  E  C".  The  following  assumptions  are  made; 


(Al)  (A,  jBa)  is  stabilizable  and  (Cj,  A)  is  detectable 

(A2)  Di2  is  full  column  rank  with  |  B12  B±  j  unitary  and  B21  is  full  row  rank  with 
unitary. 


B21 

Bi 


(A3) 

(A4) 


A  -  jul  B2 
C\  Bi2 

A  -  ju)I  Bi 
C2  B21 


has  full  column  rank  for  till  w. 


has  full  row  rank  for  all  lj. 


Assumption  (Al)  is  necessary  for  the  existence  of  stabilizing  controllers.  The  assumptions 
in  (A2)  mean  that  the  penalty  on  z  =  C\X  +  B12U  includes  a  nonsingular,  normalized  penalty 
on  the  control  u,  and  that  the  exogenous  signal  w  includes  both  plant  disturbance  and  sensor 
noise,  and  the  sensor  noise  weighting  is  normalized  and  nonsingular.  Relaocation  of  (A2)  leads 
to  singular  control  problems. 

Assumption  (A3)  relaxes  the  DGKF  assumptions  that  (Ci,  A)  is  detectable  ?.nd  ^52^1  = 
0,  and  (A4)  ’■elaxes  (A,  Bi)  stabilizable  and  BiBi^^  =  0.  Assumptions  (A3)  and  (A4)  are  made 
for  a  technical  reason;  together  with  (Al)  it  guarantees  that  the  two  Hamiltonian  matrices  in 
the  corresponding  7(2  problem  belong  to  dom(Ric).  It  is  tempting  to  suggest  that  (A3)  and 
(A4)  can  be  dropped,  but  they  are,  in  some  sense,  necessary  for  the  methods  in  this  paper  to 
be  applicable.  A  further  discussion  of  the  assumptions  and  their  possible  relaxation  will  be 
discussed  in  Section  5.2. 

It  can  be  assumed,  without  loss  of  generality,  that  7=1  since  this  is  achieved  by  the 
scalings  7~^i?ii,  7^^*C'2,  and  .  This  will  be  done  implicitly 

for  many  of  the  proofs  and  statements  of  this  paper. 


2  Preliminaries 


This  section  reviews  some  mathematical  preliminaries,  in  particular  the  computation  of  the 
various  norms  of  a  transfer  matrix  G.  Consider  the  transfer  matrix 


G(-) 


A 

B 

C 

B 

(2.1) 


with  A  stable  (i.e.,  all  eigenvalues  in  the  left  half-plane). 
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The  norm  ||G||oo  arises  in  a  number  of  ways.  Suppose  that  we  apply  an  input  w  £  C2  and 
consider  the  output  z  £  £2-  Then  a  standard  result  is  that  ||G||oo  is  the  induced  norm  of  the 
multiplication  operator  Me,  as  well  as  the  Toeplitz  operator  P+Mc  '■'H2  — >  W2. 

||(?||oo  =  sup  ||z||2  =  sup  IIP+2II2  -  sup  ||P^  Mawh 

The  rest  of  this  section  involves  additional  characterizations  of  the  norms  in  terms  of 
state-space  descriptions.  Section  2.1  collects  some  basic  material  on  the  Riccati  equation 
and  the  Riccati  operator  which  play  an  essential  role  in  the  development  of  both  theories. 
Section  2.3  reviews  some  results  on  Hankel  operators  and  introduces  the  2  X  1-block  mixed 
Hankel- Toeplitz  operator  result  that  will  play  a  key  role  in  the  Hoo  FI  problem. 

Section  2.4  includes  two  lemmas  on  characterizing  inner  transfer  functions  and  their  role 
in  certain  LFT’s  and  Section  2.5  considers  the  stabilizability  and  detectability  of  feedback 
systems. 


2.1  The  Riccati  operator 


Let  A,  iZ  be  complex  n  x  n  matrices  with  Q  and  R  Hermitian.  Define  the  2n  x  2n 
Hamiltonian  matrix 


H 


A  R 
Q  -A' 


If  we  begin  by  assuming  S  has  no  eigenvalues  on  the  imaginary  axis,  then  it  must  have  n 
eigenvalues  in  Re  s  <  0  and  n  in  Re  s  >  0.  Consider  the  two  n-dimensional  spectral  subspaces 
X^{H)  and  X+(H):  the  former  is  the  invariant  subspace  corresponding  to  eigenvalues  in  Re 
s  <  0;  the  latter,  to  eigenvalues  in  Re  s  >  0.  Finding  a  basis  for  stacking  the  basis 

vectors  up  to  form  a  matrix,  and  partitioning  the  matrix,  we  get 


X.{H)  =  Im 


X2 


(2.2) 


where  Xi,X2  €  and 


H 


Xi 

X2 


Xi 

X2 


Tx,  KeX,{Tx)<0yi 


If  Xi  is  nonsingular,  or  equivalently,  if  the  two  subspaces 


X_{H),  Im 


(2.3) 


(2.4) 


are  complementary,  we  can  set  X  :=  X2X^^.  Then  X  is  uniquely  determined  by  H,  i.e., 
H  X  is  &  function,  which  will  be  denoted  Ric\  thus,  X  =  Ric{H).  We  will  take  the  domain 
of  Ric,  denoted  dom{Ric),  to  consist  of  Hamiltonian  matrices  H  with  two  properties,  namely, 
H  has  no  eigenvalues  on  the  imaginary  axis  and  the  two  subspaces  in  (2.4)  are  complementary. 
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For  ease  of  reference,  these  will  be  called  the  stability  property  and  the  complementarity 
property,  respectively.  The  following  weU-known  results  give  some  properties  of  X  as  well  as 
verifiable  conditions  under  which  H  belongs  to  dom[Ric).  See,  for  example.  Section  7.2  in 
[Francis,  1987],  Theorem  12.2  in  [Wonham,  1985],  and  [Kucera,  1972] 

Lemma  2.1  Suppose  H  £  dom{Ric)  and  X  =  Ric(H).  Then 
fa)  X  is  Hermitian 

(b)  X  satisfies  the  algebraic  Riccati  equation 

A'X  +  XA  +  XRX  -Q  =  0 

(c)  A  +  RX  is  stable 


Lemma  2.2  Suppose  H  has  no  imaginary  eigenvalues,  R  is  either  positive  semi-definite  or 
negative  semi-definite,  and  {A,R)  is  s'abilizable.  Then  H  £  dom{Ric). 

Lemma  2.3  Suppose  H  has  the  form 

\  A  -BB'  ■ 

~  -C'C  -A' 

with  {A,B)  stabilizable  and  ranA:|  A'  +  jwl  C'  j  =  n  Vw.  Then  H  £  dom{Ric),  X  = 

Ric{H)  >  0,  and  ker(X)  C  X  :=  stable  unobservable  subspace. 


By  stable  unobservable  subspace  we  mean  the  intersection  of  the  stable  invariant  subspace 
of  A  with  the  unobservable  subspace  of  {A,C).  Note  that  if  {C,-A)  is  detectable,  then 
Ric{H)  >  0.  Also,  note  that  ker(X)  C  X  C  ker(C),  so  that  the  equation  X M  =  C  always 
has  a  solution  for  M,  for  example  the  least-squares  solution  given  by  X^C. 

We  may  extend  the  domain  of  Ric  by  relaxing  th^'  stability  requirement.  Even  if  H 
has  eigenvalues  on  the  imaginary  axis,  it  must  have  at  least  n  eigenvalues  in  Re  s  <  0. 
Suppose  that  we  now  choose  some  n-dimensional  invariant  subspace,  again  denoted  by  X_{H), 
corresponding  to  n  eigenvalues  in  Re  s  <  0  and  a  corresponding  basis  as  in  (2.2),  but  now 
satisfying 


H 


■  A'l  ■ 

X2 

Tx, 


Re  K{Tx)  <  0  V  i 


(2.5) 


This  subspace  is  not  uniquely  determined  by  H ,  but  if  it  still  satisfies  the  complementarity 
property,  then  we  can  set  AT  :=  A'aA'f  *  as  before,  if  this  X  is  also  Hermitian.  We  may 
thus  define  a  new  map  Ric,  whose  domain  dom{Ric)  wiU  be  taken  to  consist  of  Hamiltonian 
matrices  H  with  the  property  that  an  A'_(fr)  exists  satisfying  the  complementarity  condition 
and  with  the  resulting  X  :=  X2X{'^  Hermitian.  To  show  that  this  is  actually  a  map,  we  have 
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to  verify  that  X  is  uniquely  determined,  which  is  not  always  the  Ccise.  In  fact,  the  conditions 
under  which  Rtc  is  actually  a  map  are  intimately  connected  with  the  conditions  on  existence 
of  Woo  optimal  controllers.  It  turns  out  that  for  the  cases  of  interest  in  the  present  paper, 
whenever  H  is  in  dom{RT^,  the  subspace  will  be  uniquely  determined.  Thus  whenever  Ric  is 
needed,  it  will  be  a  well-defined  map,  but  this  must  be  proven.  Fortunately,  these  cases  can 
essentially  be  reduced  to  spectral  factorization  problems  and  standard  theory  can  be  applied 
(e.g.  Gohberg,  Lancaster  and  Rodman  (1986)). 

We  may  further  extend  the  domain  of  RTc  by  relaxing  the  complementarity  condition. 
The  minimal  requirement  we  will  place  on  X_{H)  is  that  (2.5)  hold  and  that 

X{X2  =  (2.6) 

is  Hermitian.  Note  that  this  condition  also  does  not  depend  on  the  particular  choice  of 
basis  taken  in  (2.2).  It  is  convenient  to  define  dom{RIc)  to  be  the  set  of  those  H  for  which  a 
subspace  X_{H)  exists  and  satisfies  (2.5)  and  (2.6).  Once  again,  the  map  Ric,  from  dom{Ric) 
to  n  dimensional  subspaces  of  (this  is  a  Grassman  manifold)  does  not  always  exist  as  the 
subspace  is  not  uniquely  determined  by  H .  _ 

The  same  remarks  about  Ric  as  a  map  apply  here  to  Ric.  These  have  been  introduced  in 
order  to  treat  the  optimal  Ctise,  but  their  use  will  be  limited  as  this  case  is  not  analysed  in 
detail.  Note  that  dom{Ric)  C  dom{Rl^  C  dom{Ric).  Also,  if  H  Q  doTn{Ric)  then  Ric  and 
Rtc  are  obviously  well-defined  maps  and  Ric{H)  =  Ric{H). 


2.2  Computing  Woo  norm 


For  the  transfer  matrix  G{s)  in  (2.1),  with  A  stable,  define  the  Hamiltonian  matrix 


H 


A  +  BR-^D'C  BR-^B' 

-C\I  -  DD')-^C  -{A+  BR-'^D’Cy 


B 

-CD 


R-^ 


D'C  B' 


(2.7) 

(2.8) 


where  R  =  I  -  D'D  .  The  following  lemma  is  essentially  from  [Anderson,  1967],  [Willems, 
1971],  and  [  Boyd  et  al.,  1989]. 


Lemma  2.4  /.  Let  a{D)  <  1,  then  the  following  conditions  are  equivalent: 

(a)  IIGIloo  <  1 

(b)  H  has  no  eigenvalues  on  the  imaginary  axis 

(c)  H  G  dom{Ric) 

(d)  H  G  dom{Ric)  and  Ric{H)  >  0  (Ric{H)  >  0  if{C,A)  is  observable) 
II.  Lei  9{D)  <  1,  then  the  following  conditions  are  equivalent: 
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(a)  IIGIloo  <  1 

(b)  H  £  dom{Rtc) 

(c)  H  6  dom{Rtc)  and  Rtc{H)  is  unique  with  Rtc{H)  >  0  (Ric{H)  >  0  if{C,A)  is  observ¬ 
able) 

Proof  From 


A 

0 

-B  1 

(7  -  G~G)(s)  =  -C'C 

-A‘  CD 

D'C 

B' 

^  J 

it  is  immediate  that  H  is  the  j4-matrix  of  (7  —  It  is  easy  to  check  using  the  PBH 

test  that  this  realization  has  no  uncontrollable  or  unobservable  modes  on  the  imaginary  ajds. 
Thus  H  has  no  eigenvalues  on  the  imaginary  axis  iff  (7  -  has  no  poles  there,  i.e., 

(7  —  €  'RCoa-  So  to  prove  the  equivalence  of  (la)  and  (Ib)  it  suffices  to  prove  that 

||G|U  <  1  o  (7  -  e  7^£cc 

I^^IIG^lloo  <  then  I -G(juj)*G(juj)  >  0,  Vu>,  and  hence  (7-6'"“<j)“*  £  TZCoo-  Conversely, 
if  ||G||oo  >  1|  then  a[G{juj)\  =  1  for  some  tu,  i.e.,  1  is  an  eigenvalue  of  G{juyG{juj),  so 
I  -  G{ju})*G{ju})  is  singular.  Thus  (la)  and  (Ib)  are  equivalent. 

The  equivalence  of  (Ib)  and  (Ic)  follows  from  Lemma  2.2,  emd  the  equivalence  of  (Ic)  and 
(Id)  follows  from  Lemma  2.1  and  standard  results  for  solutions  of  Lyapunov  equations. 

The  proof  of  part  II  is  more  involved  and  is  given  by  the  established  results  on  spectral 
factorization  as  in  Gohberg  et  a/.(1986),  since  7  -  G'^G  >  0  for  all  s  =  ju.  ■ 

In  part  II  it  was  assumed  that  a{D)  <  1  so  that  the  Hamiltonian  matrix  could  be  de¬ 
fined.  Alternatives  that  avoid  this  are  to  consider  Linear  Matrix  Inequalities  or  the  deflating 
subspaces  of  matrix  pencils.  This  is  discussed  more  in  Section  5.2.5. 

Lemma  2.4  suggests  the  following  way  to  compute  an  Hoo  norm:  select  a  positive  number 
7;  test  if  ||G||oo  <  7  by  calculating  the  eigenvalues  of  H\  increase  or  decrease  7  accordingly; 
repeat.  Thus  Tfoo  norm  computation  requires  a  search,  over  either  7  or  u.  We  should  not 

be  surprised  by  similar  characteristics  of  the  Koo-optimal  control  problem.  A  somewhat 

cLDalogous  situation  occurs  for  matrices  with  the  norms  \\M\\\  =  trace(Af*M)  and  ||M||oo  = 
a[M\.  In  principle,  HMjlj  can  be  computed  exactly  with  a  finite  number  of  operations,  as  can 
the  test  for  whether  g{M)  <  7  (e.g.  7*7  -  M*M  >  0),  but  the  value  of  a{M)  cannot.  To 
compute  a(M)  we  must  use  some  type  of  iterative  algorithm. 

2.3  Mixed  Hankel-Toeplitz  Operators 

It  will  be  useful  to  characterize  some  additional  induced  norms  of  G(s)  in  (2.1)  and  its 
associated  differential  equation 

X  =  Ax  -f  Bw 

z  =  Cx  -f-  Dw  (2.9) 
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with  A  stable.  We  will  prove  several  lemmas  that  will  be  useful  in  the  rest  of  the  paper.  It 
is  convenient  to  describe  all  the  results  in  the  frequency-domain  and  give  all  the  proofs  in 
time-domain. 

Consider  first  the  problem  of  using  an  input  w  G  £2-  to  maximize  ||P+2||^.  This  is 
exactly  the  standard  problem  of  computing  the  Hankel  norm  of  G  (i.e.,  the  induced  norm  of 
the  Hankel  operator  P+Mq  ■  —*  ^2)1  can  be  expressed  in  terms  of  the  Gramians  Lc 

and  Lo 

ALc-^  LcA' ^  BB' A' Lo  ^  L^A-^  C'C  =  Q  (2.10) 

Although  this  result  is  weU-known,  we  will  include  a  time-domain  proof  similar  in  technique 
to  the  proofs  of  the  new  results  in  this  paper. 

Lemma  2.5  sup  IIP+2II2  =  sup  \\P+MgV)\\1  -  p{LoLc) 

tue8£j_  weBH^ 


Proof  Assume  (A,  B)  is  controllable;  otherwise,  restrict  attention  to  the  controllable  sub¬ 
space.  Then  Lc  is  invertible  and  w  G  £2-  can  be  used  to  produce  any  i(0)  =  iq  given 
z(— 00)  =  0.  The  proof  is  in  two  steps.  First, 

To  show  this,  we  can  differentiate  x(t)'L~'^i(t)  along  the  solution  of  (2.9)  for  any  given  input 
w  as  follows: 

■^(x'L~^x)  =  x'L~^x  +  x'L~^x  =  x\A'L~^  -I-  L~^  A)x  ■w'B'L~^x  -h  x'L~^Bw 
Use  of  (2.10)  to  substitute  for  A'L~^  -|-  L~^A  and  completion  of  the  squares  give 

^^{x'L:^x)  =  \\wf-\\w-B'l:^xr 

Integration  from  t  =  —00  to  t  =  0  with  i(— 00)  =  0  and  x(0)  =  xq  gives 
x'qLZ^Xq  =  ||u;||^  -  lltu  -  B'L-'^x\\l  <  ||u;||i 

If  V}{t)  =  B'e~^  ^L~^xo  =  *)*io  on  (  — oo,0],  then  w  =  B'L~^x  and  equality 

is  achieved,  thus  proving  (2.11). 

Second,  given  i(0)  =  iq  and  ly  =  0,  the  norm  of  z{t)  =  Ce^^XQ  can  be  found  from 

||P+z||j  =  f  T.Qe^'^C'Ce^^xodt  =  XqLoXq 

Jo 

These  two  results  can  be  combined  as  in  Section  2  of  [Glover,  1984]; 


sup  ||P+z||^ 


sup 


WP^AW 

Mil 


XqLoXq 
=  max  /  — 

xo^O  XqIjc  Xq 


=  p[LoLc)  ■ 
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If  ||G||oo  <  1  then  by  Lemmas  2.1  and  2.4,  the  Hamiltonian  matrix  H  in  (2.8)  is  in 
dom{Ric),  X  =  Ric{H)  >  0,  >1  +  BB'X  is  stable  and 

A'X  -^XA  +  C'C  +  {XB  +  C'D)R-^{B‘X  +  D'C)  =  0  (2.12) 

Similarly,  if  a[D)  <  1  and  ||G||oo  <  1  then  by  Lemma  2.4,  the  Hamiltonian  matrix  H  in  (2.8) 
is  in  dom{Rt^,  X  =  Ric{H)  >  0,  A  +  BB'X  has  eigenvalues  in  the  closed  left  half  plane  and 
(2.12)  holds. 

The  following  lemma  offers  additional  consequence  of  bounds  on  ||G||oo-  in  fact,  this 
simple  time-domain  characterization  and  its  proof  form  the  basis  for  the  entire  development 
to  follow. 

Lemma  2.6  I.  Suppose  ||G|loo  <  1  and  i(0)  =  iq.  Then 

Slip  (ll^lli  -  ll^lli)  =  x'qXxq 

and  the  sup  is  achieved. 

II.  Suppose  that  ||G||oo  <  1,  ^{D)  <  1,  and  a:(0)  =  iq-  Then 
snp  (llzlli  -  =  x'qXxq 


Proof:  We  can  differentiate  x{t)'Xx{t)  as  above,  use  the  Riccati  equation  (2.12)  to  substitute 
for  A'X  +  X A,  and  complete  the  squares  to  get 

j^ix'Xx)  =  -||z||2  -f-  ||u;||2  -  \\R-^^^[Rw  ~  {B'X  D'C)x]\\^ 

M  w  £  C2+ ,  then  x  £  C2+ ,  so  integrating  from  t  =  0  to  t  =  00  gives 


Ml  -  ll^llz  =  _  (5'X  +  D'C)x]]\l  <  x'oXxo  (2.13) 


For  Part  I,  if  we  let  w  =  -R-^{B'X  D'C)x  =  then  w  £  C2+ 

because  A  -f  BR~^{B'X  -h  D'C)  is  stable.  Thus  the  inequality  in  (2.13)  can  be  made  an 
equality  and  the  proof  is  complete.  Note  that  the  sup  is  achieved  for  a  w  which  is  a  linear 
function  of  the  state. 

For  Part  II,  A  -f  BR~^{B'X  -t-  D'C)  may  have  imaginary  axis  eigenvalues,  hence  the 
inequality  in  (2.13)  is  still  valid,  but  may  not  give  the  supremum.  A  sequence  of  functions  Wt 
can  however  be  constructed  to  approach  the  supremum  by  considering  X^  =  Ric{H^)  where 


H,= 


B 

-CD 


{R  -I-  e^/)-' 


D'C  B' 


Then  for  in,  =  (R  +  e^I)-^B'X^  -h  D'C)x 

Ikllz  -  IKtlli  =  x'oXtXo  +  <  x'qXxo 
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Finally  taking  the  limit  as  e  — »  0  gives  the  result  by  uniqueness  oi  X  =  lime_*o  ■ 

Now  suppose  that  the  input  is  partitioned  so  that  B  =  Bi  B2  ,  D  =  Z?i  Z>2 
G{3)  =  [  Gi(s)  G2('S)  ],  and  w  is  partitioned  conformally.  Then  IIG2II00  <  1  iff 


Hw  := 


B2 

-€•02 


[  D'^C  5' 


is  in  dom{Ric),  where  R2  '■=  I  -  D'2D2-  Similarly,  0(1)2)  <  1  and  ||G||oo  <  1  iff  Hw  € 
doTn(.Ric).  In  either  case,  define  W  =  Ric{Hw),  which  will  be  unique,  and  let 

vJi  e  ,W2  e  C2^  (2.14) 

We  are  interested  in  a  test  for  sup„,ggyy  ||/+z||2  <1  (<  1),  or  equivalently 

sup  ||rtu||2  <1  (<  1)  (2.15) 

luGBW 

where  T  =  P+[Mgi  Mq,]  :  W  — *  ^^2  is  a  mixed  Hankel-Toeplitz  operator: 


10  g  VV 


■ 


-wi 

W2 


T 


Wi 

•W2 


=  P+ 


Gi  G2 


Wi 

W2 


,  W2  E  C2 


Note  that  F  is  the  sum  of  the  Hankel  operator  P+McP-  with  the  Toeplitz  operator  P+Mq^  P4.. 
The  following  lemma  generalizes  Lemma  2.4  {Bi  =  0,  Z?i  =  0)  and  Lemma  2.5  {B2  =  0,  D2  = 
0). 


Lemma  2.7  I.  (2.15)  holds  with  <  iff  the  following  two  conditions  hold: 
(i)  Hw  E  dom{Ric) 

(xi)  p{WLc)  <  1 

II.  (2.15)  holds  with  <  iff  the  following  two  conditions  hold: 

(i)  Hw  €  dom(Rt^ 

(ii)  piWLc)  <  1 


Proof  As  in  Lemma  2.5,  assume  (A,  B)  is  controUable;  otherwise,  restrict  attention  to  the 
controllable  subspace.  By  Lemma  2.4,  condition  (i)  is  necessary  for  (2.15)  for  both  cases,  so 
we  will  prove  that  given  condition  (i),  (2.15)  holds  iff  condition  (ii)  holds.  By  definition  of 
W,  if  Tu  e  W  then 

\\P^z\\\  -  Ml  =  \\PMl  -  \\P^^2\\l  -  \\PMl 

Note  that  the  last  term  only  contributes  to  ||P+z|l|  through  ®(0).  Thus  if  is  invertible, 
then  Lemma  2.6  and  (2.11)  yield 


13 


(2.16) 


sup  { l|i  +  z!ll  -  llti^lla  I  ®(0)  =  lo}  =  XqW’io  -  x'qL^  ^lo 

For  part  I  we  will  prove  the  equivalent  statement  that  p{WLc)  >  1  iff  supj^ggyv’  >  1- 

The  supremum  is  achieved  in  (2.16)  for  some  tu  6  W  that  can  be  constructed  from  the 
previous  lemmas.  Since  p(WLc)  >  1  iff  3  iq  ^  0  such  that  the  right-hand  side  of  (2.16)  is 
>  0,  we  have,  by  (2.16),  that  p{WLc)  >  1  iff  3  6  W,  u;  ^  0  such  that  ||P+2||^  >  ||x«||2-  But 

this  is  true  iff  sup„,ggyj;  lirtuHz  >  1- 

For  pcirt  II,  note  that  (2.15)  holds  with  <  iff 

sup  lirii;|l^  -  ||u;||^  <  0 
wesw 


which  by  (2.16)  is  iff  p{WLc)  <  1-  ■ 

The  FI  proof  of  Section  3.3  will  make  use  of  the  adjoint  F*  :  K2 


r*2  = 


■  P-iG^z)  ■ 

'  P.Gi  ' 

G^z 

^2 

W,  which  !.<:  givpn  by 


(2.17) 


where  P^Gz  P^(Gz)  =  {P^Mg)z.  That  the  expression  in  (2.17)  is  actuedly  the  adjoint  of 
r  is  easily  verified  from  the  definition  of  the  inner  product  on  vector- valued  £2,  expressed  in 
the  frequency-domain  as 


1  f°° 

<  xi,X2  >:=  —  J  Xi{ju)*X2{ju))dtv 


(2.18) 


The  adjoint  of  F  :  W  — »  7<2  is  the  operator  T*  :'H2  -*  W  such  that  <  2,  Fu;  >  =  <  F’2,  w  > 
for  all  tn  €  W,  2  €  7<2-  Directly  using  the  definition  in  (2.18),  we  get 

<  Z,TW  >  =  <  Z,P+{GiWi  +  G2'W2)  >  =  <  ZfGiWi  >  +  <  Z,G2‘W2  > 

=  <  P-{Gi  z),Wi  >  -\-  <  G2Z,'W2  > 

=  <  F*2,U>  > 


2.4  LFT’s  and  inner  matrices 


A  transfer  function  G  in  TZHoo,  is  called  inner  if  G~G  =  I,  and  hence  G{jw)*G{juj)  =  I 
for  all  u.  Note  that  G  inner  implies  that  G  has  at  least  as  many  rows  as  columns.  For  G 
inner,  and  any  q  e  C”*,  w  £  C2,  then  ||G(ja;)g||  =  ||9||,  Vw,  and  ||Gu;||2  =  ||in||2-  Because 
of  these  norm  preserving  properties  inner  matrices  will  be  central  to  several  of  the  proofs. 
In  this  section  we  give  a  characterization  of  inner  functions  and  some  properties  of  linear 
fractional  transformations.  First,  we  present  a  state-space  characterization  of  inner  transfer 
functions  analogous  to  Lemma  2.4  that  is  well-known  and  simple  to  verify  (see  [Anderson, 
1967],  [Wonham,  1985],  [Glover,  1984]). 


Lemma  2.8  Suppose  G  = 


■  A 

B 

C 

D 

with  {C,  A)  detectable  and  Lo  =  L'g  satisfies 
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A!U  +  loA  +  C'C  =  0. 

Then 

(a)  Lo  >  0  iff  A  is  stable 

(b)  D'C  +  B'Lo  ----  0  implies  G~G  =  D'D 

(c)  Lo  >  0,  {A,B)  controllable,  and  G~G  =  D'D  implies  D'C  +  B'Lo  =  0. 

The  next  lemma  considers  linear  fractional  transformations  with  inner  matrices  and  is 
based  on  the  work  of  Redheifer  (1960). 

Lemma  2.0  Consider  the  following  feedback  system, 


P 


Suppose  that  P-P  =  I,  PaV  G  'R’Tioo,  Q  is  a  proper  rational  matrix.  Then  the 
following  are  equivalent: 

(a)  The  system  is  internally  stable  and  well-posed,  and  |ir*u,lloo  <  1- 

(b)  Q  €  TZUoo  and  HQIU  <  1- 

Proof  (b)  =>  (a).  Internal  stability  and  well-posedness  follow  from  P,Q  £  'R-'Hoo,  ||/’22||c»  <  !> 
llQIloo  <  1,  and  a  small  gain  argument.  To  show  that  ||r*,„||oo  <  1  consider  the  closed-loop 
system  at  any  frequency  s  =  jw  with  the  signals  fixed  as  complex  constant  vectors.  Let 
liQlloo  =:  c  <  1  and  note  that  Ty,r  =  P^iil  -  P22Q)  €  RHoo-  Also  let  k  :=  HTwrlloo-  Then 
||u;||  <  «||r||,  and  P  inner  implies  that  ||z||^  +  ||r||*  =  ll^ll*  +  Therefore, 

||z|P  <  Ikll*  +  (e*  -  l)||rf  <  (1  -  (1  - 
which  implies  ||Tztu||oo  <  1  • 

(a)  =»  (b).  To  show  that  ||(?||oo  <  1  suppose  there  exist  a  (real  or  infinite)  frequency  u 
and  a  constant  nonzero  vector  r  such  that  at  s  =  jw  ,  ||Qr||  >  ||r||.  Then  setting  w  =  P2i{I- 
•^23^)1’,  V  =  Qr  gives  v  =  T^^w.  But  as  above,  P  inner  implies  that  11^11^  + llr|l^  =  |lty|P  +  llv|P 
and  hence  \\z\\^  >  ||w||*,  which  is  impossible  since  |lT,,„||oo  <  1-  It  follows  that  ^{Q{ja}))  <  1 
for  all  u,  i.c.,  ||Q||oo  <  1>  since  Q  is  rational. 

Finally,  Q  has  a  right-coprime  factorization  Q  =  NM~^  with  N,M  £  RHoo-  We  shah 
show  that  M~^  £  RHoo-  Since  To^^P^i  =  Q{I—P22Q)~^  it  has  the  right-coprime  factorization 
=  N{M  —  P22N)  '  But  since  Tov>P2i  €  RTtoo^  so  does  (M  —  ^22-^)“^.  This  implies 
that  the  winding  number  of  det(Af  -  P22N),  as  s  traverses  the  Nyquist  contour,  equals  zero. 
Furthermore,  since  det(Af  -  aP22N)  /  0  for  ah  a  in  [0,1]  and  ah  s  =  jo;  (this  uses  the  fact 
that  IIP22II00  <  1  and  IIQlloo  <  1),  we  have  that  the  winding  number  of  det  M  equals  zero  too. 
Therefore,  Q  £  RHoo  and  the  proof  is  complete.  ■ 


U/ 


1 

p  — 

■  Pll 

Pl2  ' 

V 

J  — 

P21 

P22 

£  RHc 
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2.5  LFT’s  and  stability 

In  this  section,  we  consider  the  stabilizability  and  detectability  of  feedback  systems.  The 
proofs  in  this  section  are  very  routine  and  use  standard  techniques,  principally  the  PBH  test 
for  controllability  or  observability,  so  they  will  only  be  sketched. 

Recall  the  realization  of  G  from  Section  1.4  and  suppose  that  A  6  and  that  z,  y,  w 
and  u  have  dimension  pi,  P2,  mi,  and  m2,  respectively.  Thus  Ci  G  B2  G  and 

so  on.  Now  suppose  we  apply  a  controller  K  with  stabilizable  and  detectable  realization  to 
G  to  obtsdn  .  For  the  following  lemma,  we  do  not  need  the  assumptions  from  Section  1 .4 
on  G  for  the  output  feedback  problem. 


Lemma  2.10  The  feedback  connection  of  the  realizations  for  G  and  K  is, 

=  71  +  m2  for  all  ReX  >  0. 


(a)  detectable  if  rank 

(b)  stabilizable  if  rank 


A-\I  82 
Cl  Di2 


A -XI  Bi 
C2  D21 


=  n  +  p2  for  all  ReX  >  0. 


Proof  Form  the  closed-loop  state-space  matrices  and  perform  a  PBH  test  for  controllabil¬ 
ity  and  observability.  It  is  easily  checked  that  any  unobservable  or  uncontrollable  modes 
must  occur  at  X  violating  the  above  rank  conditions  (see  Limebeer  and  Halikias  (1988)  or 
Glover(1989)  for  more  details),  hence  giving  the  results.  ■ 


3  Full  Information  and  Full  Control  Problems 


In  this  section  we  discuss  four  problems  from  which  the  output  feedback  solutions  wiU  be 
constructed  via  a  separation  argument.  These  special  problems  are  central  to  the  whole 
approach  taken  in  this  paper,  and  as  we  shall  see,  they  are  also  important  in  their  own  right. 
All  pertain  to  the  standard  block  diagram. 


but  with  different  structures  for  G.  The  problems  are  labeled 
FI.  Full  information 
FC.  Full  control 

DF.  Disturbance  feedforward  (to  be  considered  in  section  4.1) 
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OE.  Output  estimation  (to  be  considered  in  section  4.1) 

FC  and  OE  are  natural  duals  of  FI  and  DF,  respectively.  The  DF  solution  can  be  easily 
obtained  from  the  FI  solution,  as  shown  in  Section  4.1.  The  output  feedback  solutions  will 
be  constructed  out  of  the  FI  and  OE  results.  A  dual  derivation  could  use  the  FC  and  DF 
results. 

The  FI  and  FC  problems  are  not,  strictly  speaking,  special  cases  of  the  output  feedback 
problem,  as  they  do  not  satisfy  aU  of  the  assumptions.  Each  of  the  four  problems  inherits 
certain  of  the  assumptions  A1-A4  from  Section  1.4  as  appropriate.  The  terminology  and 
assumptions  will  be  discussed  in  the  subsections  for  each  problem.  In  each  of  the  four  cases, 
the  results  are  necessary  and  sufRcient  conditions  for  the  existence  of  a  controller  such  that 
ll^’zuilloo  <  7  the  family  of  all  controllers  such  that  ||T*„,||oo  <  7.  In  aU  cases,  K  must  be 
admissible. 

The  "Hoo  solution  involves  two  Hamiltonian  matrices,  Hoa  and  Jqo  which  are  defined  as 
follows: 


R  :=  Di.Di.- 
R  :=  D.iDii- 


7'/m,  0 

0  0 

7%.  o' 


where  Di,  :=  [Dn  D12] 


where  D.i  := 


^11 

£>21 


5oo 

:= 

A 

.  -qci 

0 

-A' 

- 

■  1 

B 

-CJDx. 

Joo 

;  = 

A' 

0 

~A 

- 

■  C' 
-BtDi, 

D,iB[  C 

If  Hoo  6  dom{RJ^  then  let  Xi,  X2  be  any  matrices  such  that 


■  ■ 

'  Xi  ■ 

X2 

X2 

Tx,  X[X2  =  X^Xu  ReA,(rx)<0Vz 


Similarly  if  Joo  G  dom{Rt^  then  let  Vi,  Y2  be  any  matrices  such  that 

Ty,  YlY2  =  Y^Yi  ReAi(7V)<0Vi 

Further  if  in  addition  Hoo  G  dom{RTc)  and/or  Joo  G  dom{Ric)  then  define, 
X^  :=  X2X^\  Foo  :=  FaFf* 


'  Fi  ■ 

‘  Fi  ■ 

>2 

Fa 

(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) . 


Finally  define  the  ‘state  feedback’  and  ‘output  injection’  matrices  as 


F 

L 


:=  -R-^  [Z?i.C7,  +  H'Aooj 
£2  ]  :=  +  FooC'lR-' 


(3.6) 

(3.7) 
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3.1  Problem  FI:  Full  Information 

In  the  FI  special  problem  G  has  the  following  form. 


cases,  a  suboptimal  controller  may  exist  which  uses  just  the  state  feedback  z,  but  this  will 
not  always  be  possible.  While  the  state  feedback  problem  is  more  traditional,  we  believe  that 
the  full  information  problem  is  more  fundamental  and  more  natural  than  the  state  feedback 
problem,  once  one  gets  outside  the  pure  W2  setting. 

The  assumptions  relevant  to  the  FI  problem  which  are  inherited  from  the  output  feedback 
problem  are 

(Al)  (A,  B3)  is  stabilizable. 

(A2)  D\2  is  full  column  rank  with  D12  Dx  ]  unitary. 

(A3)  has  full  column  rank  for  aU  w. 

Gi  JJ12 

The  results  for  the  Full  Information  case  are  as  follows: 

Theorem  S.l  Suppose  G  is  given  by  (3.8)  and  satisfies  Al-AS.  Then 

(a)  3K  such  that  ||r,u,||oo  <  1  €  dom{Ric),  Ric{H^)  >  0 

(b)  If  <  1  then  3K  such  that  l|r„„||oo  <  1  Ego  €  dom{Ric),  X{X2  = 

>  0-  Xi  and  X2  are  defined  in  (3.3). 

(c)  All  admissible  K{s)  such  that  ||r,„||oo  <  1  are  given  by 

forQ^UcUgg,  ligiloo  <  1. 

Note  that  the  sufficiency  proof  for  part  (b)  is  omitted.  We  will  prove  the  FI  results  and 
the  FC  results  follow  by  duality. 
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3.2  Motivation  of  the  proofs  for  Problem  FI 

We  will  first  motivate  the  proof  by  considering  a  completion  of  the  squares  assuming  that 
Xoo  >  0  and  exists.  Let  us  factor 


R  = 


r{  -I  0  1  [  Ti  0 
0  /  J  [  0  I  \[t2  I 


=:  T'JT 

^  ■  -T/Tx  +  1  ^  r  -  J 

Ti  ^  J  I  ^12^11  ^ 

=>Tj  =  DijiJii,  Tx'Ti  =  7  - 
Now  J  -  i_D'^Di\  >  0  since  at  s  =  oo 


J^i{G,K){oo)  =  73x1  +  Dx2K{oo)  ^  J  J 

and  1  >  'a{!Fi{G,  A')(oo))  >  a(73^7?ii). 

Now  consider  the  Riccati  equation  for  Xooi 

[Xoo  -l]  Hoc  =  0 

=>  X^A  +  A'X^  +  C[Ci  -  F'RF  =  0 


smd  observe  that 


(3.10) 

(3.11) 

(3.12) 


(3.13) 


=  [w'  u'j7?;,73i,  ^  -  w'w  +  x'iC[Di,  +  XooB)  ^ 

+  [  u;'  tt'  ]  (73i.Ci  +  B'X^)x  +  x'F'RFx 
=  z'z  -  x'C[Cix  -  w'w x'XooB  +  w'  u'  B'XooX  +  x'F'RFx 

U  L 

=  z'z  -  w'w  +  x'XooAx  +  x'A'XooX  +  x'XooB  ^  +  w'  u'  ]  B'XooX 

U  1 

=  z'z  -  w'w  +  -^{x'XooX) 


Integrating  from  t  =  0  to  oo  with  i(0)  =  i(oo)  =  0  gives 

lUlli  -  IIHli  =  +  .  -  [  Tj  /  ]  fx||=  -  ||r,(»  -  /■.x)||i. 


(3.14) 
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Hence  to  obtain  ||z||2  <  ||io||2  we  require  1(72^  +  “-  ^2  /  ■f’a;||2  <  ~  ■^i2:)||2 

we  see  that  in  some  sense  the  “worst  w”  is  FiX,  whereas  the  “best  u”  is  —T2W  +  T2  /  j  Fx. 
Notice  that  in  the  case  T2  ^  0  (=»  Du  7^  0)  the  natural  full  information  controller  uses  both 
w  and  X. 


3.3  Proofs  for  Problem  FI:  Necessity 

(a)  If  there  exists  an  admissible  controller  such  that  |jT,„,||oo  <  1,  then 

Hoo  G  dom{Ric),  Ric{Hca)  >  0  (3.15) 

(b)  If  there  exists  an  admissible  controller  such  that  ||Tj,j,||oo  <  1,  then 

Hoc  e  dom(AxZ).  X(A'2  =  >  0.  (316) 

We  will  prove  a  slightly  stronger  result,  but  before  that,  we  need  some  preliminary  results. 
Let  us  first  consider 

where  T  and  J  are  given  in  (3.10).  Note  that 

Di,T-^  =  [  DxD'^DuT{^\  Du 

BT-^  =  [  (Hi  -  fl2T2)T,-\  H2  ]  =:  [  Hi  H2  ]  (3.17) 

/  -  DuR-^D\,  =  /  -  H12HI2  +  Dj.D'j_DiiT{^T{-^D'uD^D'^ 

=  D^in  D'^DiriTiT,r^D[,D^)D'^ 

=  (3.18) 

where 

S  I  —  Z?(lZ?iiZ?jjZ?x  >  0. 

Hence 

N  +  HiH( 

-C[DxS-^D'^Ci  -N' 

where 

N  :=A-  BiD'uCi  +  BiT{-^ D[^D^D'^Ci 

Next  we  will  show  that  we  can  assume  without  loss  of  generality  that  the  pair  {D'^Ci,  ~N) 
is  detectable.  This  simplifies  the  technical  details  of  the  proof.  Thus  suppose  that  the  pair 
{D'^^Ci,  —N)  is  not  detectable  or  equivalently  that  {D'^Ci,  -  A  +  B2D[2Ci)  is  not  detectable. 
That  is,  [A  -  B2D12C1)  has  stable  modes  that  are  not  observable  from  D'^C\  (note  that 
modes  of  (i4  -  B2D'uCi)  on  the  imaginary  axis  are  observable  from  D'^C\  by  .43).  If  we  now 
change  state  coordinates  so  that 


~  “  -C[Cx  A' 


-C[Di. 


D'uCu  B' 
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A 

B 

L 

61. 

An 

i4i2 

6ii 

621 

= 

A21 

■^22 

612 

622 

Cn 

612 

611 

612 

with  A\2  -  B2\D'^2^i2  =  0>  =  0.  {D'^Cw,  -An  +  52il>i2^ii)  detectable  and  (>122  - 


B22D^2^'^i)  stable,  then  the  state  equations  for  the  system  with  controller  K  = 
can  be  written  as 


'  A 

6  ‘ 

C 

6 

z 

X2 

X 

U  +  -Di2^12®2 


j4iiIi  +  Biiw  +  -621(^12^12*2  +  u) 
<^11*1  +  BiiW  +  .Di2(-Dj2C'i2*2  +  ti) 
■^22*2  +  -^21*1  +  Bi2‘W  +  B22‘U 

Ax  +  BiXj  +  B2X2  +  B3W 

Cx  +  -Dill  +  62*2  +  D3W  -f  1)12^12*2 


If  the  controller,  K,  is  admissible  with  K)\\go  <  1  (<  1),  then  the  above  state 


equations  show  that  the  subsystem  Gi  = 


»ii 


'11 


B 


11 


B 
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D 


11 


D 


12 


troller,  Ki  (given  by  the  final  three  equations  above),  which  satisfies  \\J^i{G\,  K\ 
(<  1).  Furthermore,  suppose  we  can  find  a  suitable  stable  invariant  subspace 


also  has  an  admissible  con- 

00  1 


'  A-„ 
^21 


for  the  Hamiltonian  for  G\  then 

■  ^11  O' 

0  / 

X2\  0 
0  0 


will  be  suitable  for  G  since  {A22  -  6226120^12)  *s  stable.  We  will  therefore  assume  that 
{D'^C\,  -A  -f  6261261)  is  detectable  for  the  remainder  of  the  necessity  proof. 


The  proof  also  requires  a  preliminary  change  of  variables  to 
u  :=  u  -  Fqx 

This  change  of  variables  will  neither  change  internal  stability  nor  the  achievable  norm  since 
the  states  can  be  measured.  The  matrix  Fq  is  the  optimal  state  feedback  matrix  for  a  corre¬ 
sponding  'H2  problem  as  given  below.  By  Lemma  2.3  the  Hamiltonian  matrix 

Tj  A  —  626J261  —6262 

-C'6x6'^6i  -(A-B2D{2Ciy 


21 


belongs  to  dom{Ric)  since  {A,  B2)  is  stabilizable,  and  Xq  Ric[Ho)  >  0  since  -A  + 

B2D12C1)  is  detectable.  Define 


Fq  — (Di2C'i  +  B2X0),  Apo  A  +  B2FQ,  CiFi^  C\  +  D12FQ 


Gc(s) 


^Fo 

C\Fo 

Suppose  2?x  is  any  matrix  making  [D\2  D±]  an  orthogonal  matrix,  and  define 
U  U± 


■^Fo 

B2  -Xo^C{Dx  ■ 

ClFo 

^12  Dx 

Then  the  transfer  function  from  w  and  t/  to  z  becomes 
z  =  \  I  I  \=GcW  +  Ut' 


Afo 

B\  B2 

C\Fa 

Dll  Di2 

w 

V 


(3.19) 


(3.20) 


The  last  result  needed  for  the  proof  is  the  following  lemma  which  is  easily  proven  using 
t.,emma  2.8  by  obtaining  a  state-space  realization,  and  then  eliminating  uncontrollable  states 
using  a  little  algebra  involving  the  Riccati  equation  for  Xq. 


Lemma  S.2  [U  Ux\  is  square  and  inner  and  a  realization  for 


U 


is 


U  U± 


A-Fo 

B2 

-X^^C'iDx  ■ 

B'lXo  +  D'liCiFg 

DiiD\2 

D'liDx 

G  nH2 


(3.21) 


This  implies  that  U  and  Ux.  are  each  inner,  and  both  Uf^Gc  and  C/~Gc  are  in  TZH^  ■ 

We  are  now  ready  to  state  and  prove  the  main  result. 

Proposition  S.3  I.  If  sup  min  ||z||2  <  1  then  G  dom{Ric)  and  Ric[Hoo)  >  0. 

«;€B£j+  ‘'€£3+ 

II.  If  sup  min  ||2||2  <  1  then  Ego  G  dom{Rtc)  and  X{X2  =  X2X1  >  0.  Xi  and  X2 

*'€£3+ 

are  defined  in  (3.3). 


Proof  of  Proposition  Since  [ff  17x]  is  square  and  inner  by  Lemma  3.2,  ||2||2  =  ||[f/  Ux]^z\\2, 
and 


U  Ux]'^  z  = 


U'~GcU)  -f  1/ 
U^GcW 


Since  u  G  712,  its  optimal  value  is  1/  =  -P+U~GcW  and  the  hypotheses  of  the  proposition 
imply  that 
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sup 

weB7<2 


P.iU'-Gcw) 

U^GcW 


<1  (<1) 


Mixed  Hankel-Toeplitz  operators  of  this  type  were  considered  in  Section  2.3.  We  can  define 
the  adjoint  operator  F*  :  £2+  — *  VV  (W  from  (2.14))  by 


GcW 


’  P.iU'^Gcw)  ' 

P.U~  ' 

U^GcV) 

T*w  = 

of  the  operator  F  :  W  — >  W2  given  by 


92 


=  P+iG:{Uq,  +  Cfx92))  =  P+G:  [uu^ 


91 

92 


Hence 


sup  IIF9II2  <  i  (<  1) 
qeBW 

This  is  just  the  condition  (2.15),  so  from  Lemmas  2.4  and  2.7  and  (3.21)  we  have  that 

iiGrc^xiioo  <1  (<  1) 

amd  hence  Hw  €  dom{Ric)  (  or  Hw  €  dom{Iii^)  where  (substituting  for  the  Riccati  equation 
for  .X^o  and  noting  that  B{Xo  +  D'hCif^  =  T{B[Xo  +  D'nDxD'j^Ci,  see  (3.17), 


iHw)u  =  Af,  +  {-Xo^C{Dx)D'xDx,{T{Ti)-\T{B[Xo  +  D[,DxD'xC^) 
=  -Xo\A  -  H2£>'i2Ci)'Xo  -  Xo^C[DxD'xCi 
-X^^C[DxD'^DiiT-^B[Xo 
-Xo^C[DxS-'^D'^DrrD[,DxD'^C^ 

=  -Xq^N'Xo  -  Xo^C[DxS-^D'^Ci 

{Hw)i2  =  Xo^C[DxS-'^Ij\^\Xo^ 


{Hw)2i  =  -{XoBiT^+C[DxD'xDn){T[T^)-\T[B[Xo-^D[^DxD'xC^) 
=  -XqB^B[Xq  -  N'Xq  -  XoN  +  XoB2B'2Xo  -  C{DxS-^ D'^Ci 

It  is  now  immediate  that 

Hw  =  T-^HooT  where  T  = 


-I 

v-i 

-^0 

T’-l  _ 

0  -Xq-'  ■ 

-Xo 

0 

,  j.  — 

Xo  -I 

The  appropriate  stable  invariant  subspace  for  Hw  will  be  Im 
Hod  will  be 


I 

W 


and  hence  that  for 
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I 


I'^J’  [ 

Moreover  Lemma  2.7  will  give  that  p{WXq^)  <  1  (<  1)  and  hence  Xq  >  W  {Xq  >  W) 
giving  that 

(J  -  Xo^WyXo  =  Xo  -  >  0  (>  0) 


Xoo  =  Xo{Xq  -  wy^Xo  >  0 

in  case  (a).  This  completes  the  necessity  proof  for  both  parts  (a)  and  (b). 


3.4  Proofs  for  Problem  FI:  Sufficiency 
All  admissible  X(s)  such  that  ||Ti,u||oo  <  1  are  given  by 


for  Q  e  UcHoo,  ligiloo  <  1. 

Note  that  this  contains  the  t/part  of  (a). 

Before  beginning  the  proof,  we  will  perform  a  change  of  variables  suggested  by  Section  3.2. 
Firstly  change  the  input  variable  to 

V  =  u  +  T2'w  -  I  72  f  Fx 

with  the  corresponding  controller 

iftmp(^)  =  K{S)  +  [  -  [  T2  /  ]  F  T2  ] 

and  state  equations 

X  =  Apx  +  {Bi  -  B2T2)w  +  B2V 
z  =  Cipx  +  D  ^D'^D\\vi  +  D\2V 

where 

Ap  :=  [A  +  B2  T2  I  F);  C\f  =  C7i  +  D\2  T2  /  J  F 

Also  define  the  new  feedback  variable 
r  ;=  71(1^  -  Fii) 

Now  suppose 
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=  Qis)Tr  [  -Fi  /  ]  ; 


that  is 

V  =  Qr 

This  gives  the  following  feedback  configuration  in  which  one  would  expect  from  (3.14)  that 
=  I  since  —  H^Hi  =  ll^lla  —  IIt’Hq  and  this  is  now  proven  together  with  the  stability 
of  Af. 


Af 

B\  —  B2T2 

B2  1 

CiF 

D±D'^Dii 

D\2 

.  -TiFi 

Ti 

0 

Lemma  8.4  P  £  Tlc'Hoo  ,  P'"P  —  I  and  €  TZcHoo- 
Proof  The  observability  Gramian  of  P  is  Xoa  since 

A>pX^  +  X^Af  +  C\fCif  +  F'  [  Ti  0  ]  F 

=  A'Xoo  +  XooA  +  C'iCi  +  F^  ^2  {b’2X^  + 

Hx^B2  +  c{d,2){t2  i]f 

1]] 

= 


(3.22) 


where  we  have  used  the  identity  -Bl^X^o  -  =  [^2  ^  \  F-  Furthermore,  since  X^o  >  0 

and  (Fi,  Af)  is  detectable  (ncte  i4F+(Fi-F2r2)Fi  =  A-\-BF  is  stable  since  X^o  =  Ric(iy«,)) 
we  have  that  Af  is  stable  by  Lemma  2.8(a).  Also 

■  DiiDxU'i  Ti'  ]  [  +  ^12  [  F2  /  ]  F  1  [  bJ  -  TiB'2  „ 

0  J  [  -  [  Ti,  0  ]  F  FJ  J  ~ 


0  I 


(£>i.(7i  +  FF+F'A-„) 
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Hence  by  Lemma  2.8(b), 


T{ 

D’12  0 

/  0  ’ 

0  I 


DxD*^Dii  Di2 
Ti  0 


as  claimed.  It  is  also  easily  shown  that  6  'J^'Hoo  since  its  poles  are  Ai(yl  -4-  BF).  ■ 

The  proof  of  sufficiency  for  Theorem  3.1(a)  and  the  class  of  all  controllers  given  in  Theo¬ 
rem  3.1(c)  can  now  be  completed.  Let  K  be  any  admissible  controller  such  that  ||rs,„||oo  <  1- 
Then  €  HHoo  and  =  Pn  -f-  Now  define  Q  =  {I  +  T^P^i  P22)~^'F^P2i  so 

that  Q{I  —  P22Q)~^P2\  =  Tvw  and  T^ui  =  Pii  +  PiiQil  -  P22Q)~^ P^i-  Since  P22  is  strictly 
proper  all  the  above  are  weU-posed  ind  Q  is  real-rational  and  proper.  Hence  Lemma  2.9 
implies  that  Q  €  TlHoo  with  !!Q||oo  <  1-  This  verifies  that  all  transfer  functions  and 
hence  T„,„,  can  be  represented  in  this  way.  ■ 


Remark- 


In  the  optimal  case  of  part  (6)  the  proof  of  sufficiency  is  more  delicate  and  to  iUustrate  the 
difficulty  the  following  example  is  given.  Let 


G  = 


1 

1 

1 

“ 

1 

0 

0 

0 

0 

1 

1 

0 

’  0  ' 

. 

0 

1 

0 

then, 


1 

0 

■  ' 

'  0  ' 

-1 

-1 

> 

X2 

1 

Tx  =  -1,  X[X2  =  0  >  0. 


An  optimal  controller  is  given  by 

u  =  Fx  -  w,=^  X  =  (F  1)1,  zi  —  X,  zz  =  V.  =  Fx  -  w, 
where  F  -f-  1  <  0  but  F  is  otherwise  arbitrary.  Clearly  for  this  controller  i  =  0  and  hence 

2i  —  0,  Z2  =  —vj. 

If  the  controller  for  the  suboptimal  case  with  7"^  =  1  -  is  applied  (see  DGKF  item 
FI.5),  then. 


Xr^  = 


i-hv/r+^ 


K{s)  =  [  -  g(s)(l  -  e'^)X^  Q{s)  ] 

An  admissible  optimal  controller  is  obtained  as  t  -♦  0  iff  Q(s)  =  -1,  in  which  case  X(s) 

'  -(1  +  ^/TT^)  -1  ' 
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3.5  Problem  FC:  Full  Control 
The  FC  problem  has  G  given  by, 


i 

A 

Bx 

10 

1 

G{s)  = 

Cl 

Dxx 

0  / 

C2 

D21 

0  0 

(3.23) 


and  is  the  dual  of  the  Full  Information  case:  the  G  for  the  FC  problem  has  the  same  form  as 
the  transpose  of  G  for  the  FI  problem.  The  term  Full  Control  is  used  because  the  controller 
has  full  access  to  both  the  state  through  output  injection  and  to  the  output  z.  The  oidy 
restriction  on  the  controller  is  that  it  must  work  with  the  measurement  y.  The  assumptions 
that  the  FC  problem  inherits  from  the  output  feedback  problem  are  just  the  dual  of  those  in 
the  FI  problem: 


(Al)  (C'2,A)  is  detectable 
(A2)  D21  is  full  row  rank  with 


D21 

D± 


unitary. 


(A4) 


A  -  ju)I  B\ 
C2  D21 


has  full  row  rank  for  all  a>. 


Necessary  and  sufHcient  conditions  for  the  FC  case  are  given  in  the  following  corollary. 
The  fartuly  of  all  controllers  can  be  obtained  from  the  dual  of  Theorem  3.1  but  these  will  not 
be  required  in  the  sequel  and  are  hence  omitted. 


Corollary  3.5  Suppose  G  is  given  by  (3.23)  and  satisfies  Al,  A2  and  A4-  Then 

(a)  3K  such  that  ||Tit„||oo  <1  Joo  ^  dom{Ric),  Ric{Joo)  >  0 

(b)  3K  such  that  ||r„i,||oo  <1  «foo  €  dom{Rtc),  X{X2  =  X2X1  >  0.  Xi  and  X2  are 
defined  in  (3.4). 


4  Main  Results:  Output  feedback 

The  solution  to  the  Full  Information  problem  of  section  3  is  used  in  this  section  to  solve  the 
output  feedback  problem.  Firstly  in  Theorem  4.2  a  so-called  disturbance  feedforward  problem 
is  solved.  In  this  problem  one  component  of  the  disturbance,  ■W2,  can  be  estimated  exactly 
from  y  using  an  observer,  and  the  other  component  of  the  disturbance,  wi ,  does  not  affect  the 
state  or  the  output.  The  conditions  for  the  existence  of  a  controller  satisfying  a  closed-loop 
Tfoo-Qorm  constraint  is  then  identical  to  the  FI  case. 

The  solution  to  the  general  output  feedback  problem  can  then  be  derived  from  the  trans¬ 
pose  of  Theorem  4.1  (CoroUairy  4.3)  by  a  suitable  change  of  variables  which  is  based  on  X^o 
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and  the  completion  of  the  squares  argument  given  in  Section  3.2  and  the  characterization  of 
all  solutions  given  in  Section  3.4. 

The  main  result  is  now  stated  in  terms  of  the  matrices  defined  in  section  3  involving  the 
solutions  of  the  Xoo  and  Voo  Riccati  equations  together  with  the  “state  feedback”  and  “output 
iujecticn”  matrices  F  and  L.  It  will  further  be  convenient  to  additionally  assume  unitary 
changes  of  coordinates  on  w  and  z  have  been  carried  out  to  give  the  following  partitions  of 
D,  Fi  and  Li. 


Theorem  4.1  Suppose  G  satisfies  the  assumptions  AI-A4  of  section  1.4- 

(a)  There  exists  an  admissible  controller  K(s)  such  that  \\Fi{G,  A')||oo  <  7  (i.e.  ||r,u-||oo  < 
7^  »/  and  only  if 

(i)  7  > 

(ii)  Boo  G  dom(Ric)  with  Xgo  =  Ric{Hoo)  >  0 
(Hi)  Joo  €  dom{Ric)  with  Yoo  =  Ric{Joo)  >  0 
(iv)  p(X^Y^)  <  72. 

(b)  Given  that  the  conditions  of  part  (a)  are  satisfied,  then  all  rational  internally  stabilizing 
controllers  K{s)  satisfying  ||.?i(Gr,ir)||oo  < 'y  are  given  by 

K  =  Fi{Ka,^)  for  arbitrary  $  £  BcHoo  such  that  ||#||oo  <  7 

where 


■  A 

Bi 

62 

Ka  =  Cl 

Dll 

Di2 

.C2 

D21 

0 

Dll  =  -Dii2iD[m{’y^I  -  DiiiiD[m)~^ D1112  -  D1122, 

Di2  G  C”**’'”**  and  D21  G  are  any  matrices  (e.g.  Cholesky  factors)  satisfying 

D12D12  —  J  -  Dii2i{,'y^I  -  D[iiiDiiii)~^ D[i2ii 
D21D21  =  /  - 

and 

62  =  2^(82  +  Li2)Di2, 
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Cj  =  -^2l(C'2  +  i^l2)) 

Bi  =  -Z^L2  +  B2D12 
Cl  =  F2  +  DuD^iCz, 
i  =  A  +  5F  +  Bibi^C2, 


where 


Zoo=(/-7‘’noA'oo). 


(Note  that  if  Du  =  0  then  the  formulae  are  considerably  simplihed.) 

The  proof  of  this  main  result  is  via  some  special  problems  that  are  simpler  special  cases 
of  the  general  problem  and  can  be  derived  from  the  FI  and  FC  problems.  A  separation  type 
argument  can  then  give  the  solution  to  the  general  problem  from  these  special  problems.  It 
can  be  assumed,  without  loss  of  generality,  that  7  =  1  since  this  is  achieved  by  the  scalings 
7~^Dii,  7'^/^C'i,  7^/^J?2,  7^^*C2,  'y'^Xoo,  7~^Foo  and  7~^iif.  All  the  proofs  will 

be  given  for  the  case  7  =  1. 

4.1  Disturbance  Feedforward 

In  the  Disturbance  Feedforward  problem  one  component  of  the  disturbance,  wi,  does  not 
affect  the  state  or  the  output.  The  other  component  of  the  disturbance,  W2  (  and  hence 
the  state  i),  can  be  estimated  exactly  from  y  using  an  observer.  The  conditions  for  the 
existence  of  a  controller  satisfying  a  closed-loop  ?<oo-norm  constraint  is  then  identical  to  the 
Full  Information  case. 

Theorem  4.2  (Disturbance  Feedforward) 

Theorem  4.I  is  true  under  the  additional  assumptions  that 

Bib'^  —  0;  A  —  B1D21C2  is  stable.  (4.2) 

In  this  case, 

Too  =  0,  Z  =  /,  i  =  -  [  0,  BiD'2i  ] 
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Proof 

(a)  The  necessity  of  the  conditions  is  immediate  from  Theorem  3.1  since  the  existence  of  an 
output  feedback  controller  implies  the  existence  of  a  state  feedback  controller.  Further,  the 
additional  condition  <  1  is  clearly  necessary  by  considering  s  =  oo.  Theorem  3.1 

also  shows  that  all  controllers  satisfying  ||.?i(G,  Jir)jjoo  <  1  are  given  by 

u  =  Q{s)Ti{w  -  Fix)  +  T^^Fix  -  tu)  +  F2X 
r  =  Ti(tu  -  Fix) 

V  =  Qr 

For  any  Q  G  RHoo,  ||<?||oo  <  1-  Also  the  transfer  function  T^w  is  obtained  from  the  block 
diagram 


and  hence 

ti  =  ((/  —  QP22)~^QP2i  ~  ^2)“^  +  T^r^x  +  F2X  (4.4) 

We  need  to  find  a  Q(s)  that  can  be  written  as  an  output  feedback.  The  assumption  of 
(4.1)  will  give  the  following  realization  for  G, 


Hence  wi  affects  z  but  neither  x  nor  y  and  we  must  firstly  find  a  (?(s)  in  (4.4)  such  that 
Tuwi  is  zero.  Since  is  zero  we  need 

=  (/  -  QP22)-\Q{P2i  +  P22T2)  -  T2)b\  =  0  (4.5) 

Using  the  state  space  realization  that  for  P21  P22  in  (3.22)  gives 

[P21  +  P22T2]b'^  =  Tib'^ 

=»  QTib'^  =  Tib'^  (4.6) 

Again  without  loss  of  generality  we  can  assume  that 
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where 


Tn  T\2 

0  Ti3 


T[T^  =  I-D[^DxD'^Dir 

and  hence 

^n^ii  =  ^  -  ^ii  11-^^1111 
T^D'x  —  D[2DiiD'^^  =  Dll  21 
T{iTi2  =  -DiiuDiii2 

^13^13  =  -  Diji2DiiiiTj^^rjj^DiiiiDiii2  -  ^1112^1112  =:  £>21^21 

Hence  (4.6)  implies  that  for  <3  =  J  Qi  <?2  |i 

Qi  =  Dii2i3’iV 

and  QQ~  <  I  implies  that 

Q2Q2  ^  ~  Dii2i(rj[iTii)  ^Dji21 

=  +  Dll 2i(/ -  DiijjDii  11  -  Dji 21-^11 21)  ^•f^li2i)  ^ 

■•=  D12D12 

where  the  indicated  inverses  exist  by  (aj(ij.  Hence 
Q2  =  D12Q3  for  Qz  G  RHooi  IIQ3II00  <  1- 

We  have  hence  shown  that  all  controllers  can  be  written  as  feedback  from  'W2  and  x  by 
substituting  for  Q  into  (4.4)  as 


+  D1121,  D1122  {FiX  —  w)  +  F2X 


=  DizQzDzii'Wz  —  DziFix)  +  Dii2ilij^[ri2u;2  —  [7’iiri2]Fii] 

+  [  D1121,  D1122  ]  Fix  -  D1122W2  +  F2X 

=  Di2Q3D2i(W2  -  D2lfll) 

+  (-Dii2i(/  -  DiiiiDiiii)~^DjiiiDiii2  -  Dii22)('U'2  -  D21D1I) 

+  ■f2a! 

=  (Dll  +  Di2Q3D2l)(n'2  -  ^12®)  +  ^2®. 
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This  gives  the  complete  family  of  controllers  in  terms  of  i  and  W2.  The  disturbance,  W2, 
and  state  z  can  be  exactly  estimated  from  the  measurement,  y,  by  means  of  an  observer  as 
follows, 

X  =  Ax  +  Bi2‘W2  +  B2U 

W2  =  -€2^  +  y 

u  =  -  i^iaz)  +  D\2V  +  ^2® 

p  =  Qzq 

q  —  D2i‘W2  -  D2\Fii^ 

It  follows  that 

(z  -  z)  =  (.4  -  5i2C'2)(z  -  z) 

and  hence  for  z(0)  =  z(0)  =  0,  z(t)  =  z(t)  and  w{t)  =  w{t)  for  all  t  >  0.  Furthermore, 
internal  stability  will  follow  from  the  stability  of  A  —  B\2C2- 

Finally  it  is  straightforward  to  verify  that  this  family  of  controllers  corresponds  exactly 
to  those  of  Theorem  4.1  with  Yoo  =  0,  Z  =  I,  and  since 

L  =  -BiD'.^R-^  =  -  [  0  512 

and  0  =  Ric(Joo)- 


The  transpose  of  Theorem  4.2  can  now  be  stated  to  obtain  another  special  case  of  Theo¬ 
rem  4.1. 


Corollary  4.3  (Output  Estimation) 

Theorem  4-i  is  true  under  the  additional  assumptions  that 

D'^Ci  =  0,  A  —  is  stable. 


In  this  case 


Xoo  =  0,  z  =  l,  F=- 


"  1 


4.2  Converting  Output  Feedback  to  Output  Estimation 

The  output  feedback  case  when  the  disturbance,  w,  cannot  be  estimated  from  the  output  is 
reduced  to  the  case  of  Corollary  4.3  by  a  suitable  change  of  variables.  Since  we  showed  in 
(3.14)  that 
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where 


V  =  u  +  T^w  —  T2,  I  Fx 

r  =  Ti(t«  -  F\x) 

We  will  perform  the  change  of  variables  with  v  replacing  z  and  r  replacing  w.  Hence 
X  —  (A  +  B\F\)x  +  B\T-y  +  H2U 

V  =  u  +  T2Tf^r  -  F2X 

y  =  C2X  +  D2\T^^r  +  D21F1X 


and  the  transfer  matrix  from 


is 


Similarly  substituting  v  for  u  in  the  equation  for  G  gives  that  the  transfer  function  from 
^  ^  ^  to  ^  ^  ^  is  P  as  defined  in  (3.22).  We  can  show  with  a  little  algebra  the  equivalence 
of  the  first  two  of  the  following  block  diagrams,  with  T^r  given  by  the  third  one. 


Lemma  4.4  Let  G  satisfy  AI-A4,  and  assume  that  Xoo  exists  and  Xoo  >  0-  Then  the  fol¬ 
lowing  are  equivalent: 

(a)  K  internally  stabilizes  G  and  ||P/(G,  A')||oo  <  1, 

(b)  K  internally  stabilizes  Gvyru  ||P<(G,;yru>  ^k^)||oo  <  1, 

(c)  K  internally  stabilizes  Gtmp  and  ||Pr(Gtnip,  .f<^)||oo  <  1, 
where  G^yru  w  given  by  (4-7)  and 


B2 

^tmp  •  — 

-D12F2 

Dn 

Di2 

C2  4-  D21F1 

D21 

0 

Proof 


(a)  ^  (b)  Referring  to  the  above  block  diagram  for  P  and  it  is  seen  by  Lemma  2.9 
that  €  RHoo  with  ||r„„||oo  <  1  iff  T„r  6  RHoo  with  ||r„r|loo  <  1-  (RecaU  that  =  I, 
P  €  RHoo,  and  G  RHoo)-  In  order  to  prove  internal  stability  of  both  systems  we  note  that 
this  is  equivalent  to  the  realizations  being  stabilizable  and  detectable.  The  realization  of  T„  is 
detectable  since  the  system  zeros  of  ((3vyru)i2  are  the  eigenvalues  of  A-\-BF  (see  Lemma  2.10). 
Further  the  realisation  of  T^r  is  stabilizable  from  r  iff  the  realisation  of  Tzw  is  stabilizable  from 
w  since  they  are  related  by  state  feedback.  Finally  if  the  realisation  of  T^r  is  internally  stable 
with  ||T„||oo  <  1  then  the  above  block  diagrams  for  =  P t{P,  Pi{G vyru,  R))  is  internally 
stable  by  a  smadl  gain  argument  and  hence  so  is  that  for  P'i{G,  K). 

(b)  (c)  Internal  stability  of  both  systems  is  equivalent  since  the  closed-loop  i4-matrices 

are  identical.  Further  note  that 

_  _  Di2  0  0  Ti  0 

i^tmp  -  1^  Q  Q  ^  J  1^  J  [  0  ^ 

cind  recall  that 

TiTi  =  I-D[,D^D'^Dn. 

Hence 

Ft{Gtmp,K)=  [  Du] 

and 

/  -  K)~Pt{Gizop,  K)  =  T[{I  -  Ti{Go^,z,  K)-Pi{Goyru,  K))T, 

hence  giving  the  equivalence  of  (6)  and  (c). 


The  importance  of  the  above  constructions  for  G^yru  and  Gtmp  is  that  they  satisfy  the 
assumptions  for  the  output  estimation  problem  (CoroUary  4.3)  since  A  -f  BF  is  stable.  Hence 
we  are  now  able  to  prove  Theorem  4.1. 


Proof  of  Theorem  4.1  (Output  Feedback) 

(a).  The  necessity  of  the  conditions  will  be  first  proved.  Let  /f  be  a  proper  controller  satisfying 
||T»u)||oo  <  1.  then  the  controller  K  |  C2  D21  j  solves  the  full  information  problem  and  hence 

D 

(it)  holds.  Similarly  K  solves  the  full  control  problem  and  hence  [Hi)  holds.  From 

D\2  J 

Lemma  4.4  K  stabilizes  Gtmp  with  llPi(Gtmpi  R^)!!*  <  li  which  satisfies  the  assumptions  for 
the  output  estimation  problem  of  Corollary  4.3,  since  AA  implies  that 


rank 


A  +  B\Fi  —  jLjI  Bi 

C2  +  DiiFi  D21 


=  n  +  mi 


Hence  we  require 

•^tmp  ^  d07Tl(^ Ric'^ 


and 


Vtmp  • —  B.ic(^JtTap)  ^  0 


where 


'tmp  — 


We  claim  that 


Ja- 


i 

0 


A'  +  F[B[  0 
-fiiH'  -A-ByFi 


■  -FiD[, 

C'2  +  FiD'2,  ■ 

R-^ 

BuB'i  —D12F2 

f?21^1  C2  +  D2\F\ 

Joo 


/  ^00 

0  I 


'trop 


where  Joo  was  defined  in  (3.2)  as 


Joo  := 

A' 

-BiB{ 

0 

-A 

- 

C 

-B,D',, 

k-^ 

D,^B[  C 

To  verify  this  claim  let 


M 

E 

N 


-D12F2 
C2  +  D21F1 

D,iB'iX„  +  C-M 
Bi,F  +  Cl 


Substituting  for  fii^oo  from  (3.7)  gives 


B[X^  = 

F,  -  D[,N 

E  = 

D,i{Fi- D[,N)  + 

Cl  +  D12F2 
-D21F1 

N 

[  -£>2ijPn 

and  hence 

(4.8) 


(4.9) 


(4.10) 
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-  R-'^E  = 


N 


(4.11) 


Now  consider  the  claim  component  by  component.  Clearly 

(<>^0)21  =  {Joo)2\  -  (<^tmp)21 

Secondly 

(Ja)22  =  -A- ^  E) 

=  -A  + BiD',^R-^M  -  Bx{B[Xc^  + D'^^N) 

—  (‘l^tmp)22 

by  (4.10). 

Finally 

(Jtx)l2  —  (*^tmp)l2  =  X03A  +  A'Xoo  +  XooB\  B[Xoo 

-  {M'  +  E*)R-\M  +  £;)  +  M'R-^M 

Substitute  from  (3.13): 

XooA  +  A'Xo^  =  -CjCi  +  F'RF 

=  -{N'  -  F'D[,){N  -  DuF)  +  F'D[,D^,F  -  F{F, 

=  -N'N  +  F'D[,N  +  N'DuF  -  F{Fi 

Equation  (4.10)  gives 

XcoBiB[X^  =  (FI  -  N'BuXFi  -  D[,N) 

and  (4.11)  and  (4.8)  give 

-M'R-^E  -  E'R-^M  -  E'R-^E  =  -  N' DuF2 

-  -  I)N 

Adding  these  three  expressions  gives  (Ja)i2  =  (•^tmp)i2  and  the  claim  that  Ja  =  -/imp  is 
verified. 

Since 

Joo  y  =  y  (AI^C'VX 

*00  ^00 

we  have 

T  f  ~  -^00^00  _  ^  ~  XooYoO  f  a'  ^  r'\ 

Jtmp  y  -  y  (A  +  L,  ) 

•*00  ^00 

and 
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^mp  :=  Ric(  Jtmp)  =  Y„{I  -  X^Y^)-^  >  0. 

It  is  readily  verified  that  this  implies  and  is  implied  by  (iv),  that  p{XooYoo)  <  1.  To  see 

y  Q 

this,  consider  Yoo  =  0  ’  ^  note  that  —  Xgoi  >  0;  conversely  note 

that  XooYoo  =  (-^^  +  YtjapXoo)~^YtrapXoo  and  hence  Timp  >  0  implies  p{XooYoo)  <  1- 

Therefore  the  necessity  of  the  condition  is  proven.  Sufficiency  also  follows  immediately 
because  of  the  equivalence  of  the  G  and  <?tmp  problems. 


(b)  Characterization  of  all  solutions 


To  characterize  all  controllers  for  G  we  just  need  to  characterize  all  controllers  for  Gtmp  using 
Corollary  4.3,  with  Timp  =  YoaZ'^^  where 

^oo  :=  (I-YooX^) 

Limp  =  -{BiD',^+YimpM')R-'^ 

=  +  Y^{-X^B^D',^  +  M'))R-^ 

=  +  Y^{C'  -  E'))R-^ 


ff'  0 


!  -^tmp  —  Oi  ^tmp  —  L 


We  can  nov.  substitute  in  the  formulae  of  Theorem  4.1  with  Gtmp  and  the  above  (•)tmp 
values  to  obtain  the  class  of  controllers. 


62  —  {B2  Z^  L\Di2  -  Z^YooN'D\2)D\2 

=  Z^\B2  -  Y^X„B2  +  LiZ)i2  -  Yoo(F'Z?;.  +  C[)D^2)bx2 

—  ^00  {B2  +  L\Di2)bi2 

by  (3.7).  The  expressions  for  C\,  62,  B\  and  A  are  then  obtained  by  a  direct  transcription 
of  the  above  expressions  and  are  hence  omitted.  This  completes  the  proof.  ■ 


5  Generalizations 

In  this  section  we  indicate  how  the  results  of  section  4  can  be  extended  to  more  general  cases. 
Firstly  the  optimal  case  is  considered  when  a  variety  of  new  phenomena  are  encountered. 
Secondly  the  removal  of  assumptions  A1-A4  is  discussed.  Finally  some  comments  are  included 
for  the  case  when  the  optimal  Tfoo-norm  is  necessarily  achieved  at  s  =  00. 

5.1  The  Optimal  Case 

In  the  optimal  case  any  combination  cf  the  conditions  of  Theorem  4.1(a)  may  be  violated.  In 
order  that  the  Hamiltonian  matrices  Hoo  and  Joo  can  be  defined  we  will  assume  that  condition 
(a)(i)  in  Theorem  4.1  is  satisfied  and  will  state  the  result  proven  in  Glover  et  al.  (1989). 
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Firstly  if  Hoo,  Joo  G  dom(Rtc)  then  there  exist 


X2 


satisfying  equation  (3.3)  and 


Y  J  satisfying  equation  (3.4).  In  the  optimal  case  Xi  and/or  Yi  may  be  singular  so 

that  Xoo  :=  and  1^00  :=  may  not  exist,  and  if  these  inverses  exist  ;= 

/  _  ‘y~^YooXoo  may  be  singular.  In  order  to  avoid  taking  these  inverses  we  will  modify  the 
definitions  of  the  ‘state-feedback’  matrix,  F  in  (3.7),  and  the  ‘output  injection’  matrix,  L  in 
(3.7),  as  follows. 


.F°  :=  -R-^[D[,CiXi  +  B'Xi] 
L°  ;=  -[Y{BiD',^+Y^C]R-^ 


Furthermore  as  in  (4.1)  we  assume  that  D,  F°,  and  L°  have  been  transformed  and  partitioned 
as  follows. 


po  1 

po  1 
^12 

F2  1 

pot 

TO  t 
^11 

Dim 

^1112 

0 

D 

D1121 

D1122 

I 

[  LI' 

0 

I 

0 

The  solution  to  the  output  feedback  problem  in  the  optimal  case  can  now  be  stated  (Glover 
et  ai  (1989)). 


Theorem  5.1  Suppose  G  satisfies  the  assumptions  AI-A4  of  section  I.4  and 

7  >  max(CT[Z?nii.  £>1112,], ^[^'1111. 

(a)  There  exists  an  admissible  controller  K {s)  sxich  that  \\Ft{G,K)\\oo  <  7  (i-e.  ||r*u;||oo  < 
7/1/  and  only  if 


(i)  Hfx>  G  dom{Rt^  with  Xi,  X2  satisfying  (3.3)  such  that  X[X2  >  0. 

(ii)  Joo  G  dom{Rtc)  with  Yi,  Y2  satisfying  (3-4)  such  that  Y(Y2  >  0. 


(iH) 


X!,X^  -I-^X'2Y2 

'y-^YfX2  YfY^ 


>  0 


(b)  Given  that  the  conditions  of  part  (a)  are  satisfied,  then  all  rational  internally  stabilizing 
controllers  K{s)  satisfying  ||.^/(<?,  ir)|loo  £  7  given  by 


K  =  Fi{Ka,^)  for  arbitrary  $  G  'R-c'^Co 


such  that  ||$||oo  <  7.  det(7  -  (iira)22(<»)#(oo))  ^  0. 


where 


Ko  := 


b\\  bi2 
b2\  0 


+ 


Cf 

C°2 


[sE-A)*\Bt  B 
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*  denotes  a  suitable  pseudo  inverse,  Dij  are  defined  in  Theorem  4.1-  and 

:=  {Y^B2  +  Ll^)Di2 
:=  -D2iiC2Xi  +  F?^) 

B\  :=  — i/5  + 

Cl  :=  Fl  +  bxib2^Cl 
A°  :=  ETx-^Btb:i^Cl  =  Tl^E^BlbliCl 
E  := 

■ 

The  descriptor  form  of  the  equations  for  the  controllers  has  been  used  as  proposed  for 
optimal  Hankel-norm  approximation  by  Safonov  et  at.  (1987).  At  optimality  E  will  typically 
be  singrdar  and  the  state-space  equations  of  Theorem  4.1  are  not  possible.  Moreover  the 
matrix  {sE  —  A°)  may  be  singular  for  all  s,  but  the  transfer  function  Ka{s)  nevertheless 
remains  uniquely  defined.  The  condition  that  det(/  —  (iira)22(oo)#(oo))  ^  0  is  required  so 
that  this  LFT  is  well-posed.  It  is  often  the  case  that  all  the  controllers  can  be  characterized  by 
$  =  M\i\M2  for  non-square  constant  matrices,  M{Mi  =  I  and  M2M2  =  I,  with  #1  G  'R-c'Hoo 
such  that  11^1  ||oo  <  7- 

The  optimal  case  may  also  occur  when  Hoo  or  Joo  have  eigen- values  on  the  imaginary  axis 
but  Hoo,  Joo  €  dom{Rtc).  In  this  case  Theorem  5.1  can  give  regular  state-space  equations 
with  E,  Xi,  and  Yi  all  invertible.  The  stable  invariant  subspace  of  Hoo  or  Joo  will  only  be 
unique  when  the  additional  constraint  that  X{X2  and  Y1Y2  are  Hermitian  is  included,  and 
this  requires  some  special  purpose  algorithms  (see  Section  5.2.5). 

5.2  Relaxing  Assumptions  A1-A4 

5.2.1  Relsocing  AS  and  A4 

Suppose  that, 

■  0  0  1  ■ 

G=  0  0  1 
110 

which  violates  both  A3  and  A4  and  corresponds  to  the  robust  stabilization  of  an  integrator. 
If  the  controller  u  =  —  ci,  for  c  >  0  is  used  then 

Ttui  —  T  »  with  ||Txu,||oo  =  f 
s  -r  c 

Hence  the  norm  can  be  made  arbitrarily  small  as  c  — »  0,  but  e  =  0  is  not  admissible  since  it 
is  not  stabilizing.  This  may  be  thought  of  as  a  case  where  the  7foo-optimum  is  not  achieved 
on  the  set  of  admissible  controllers.  Of  course,  for  this  system,  Hoo  optimal  control  is  a  silly 
problem,  although  the  subop timal  case  is  not  obviously  so. 


If  one  simply  drops  the  requirement  that  controllers  be  admissible  and  removes  assump¬ 
tions  A3  and  A4,  then  the  formulae  in  this  paper  will  yield  u  =  0  for  both  the  optimal 
controller  and  the  suboptimal  controller  with  $  =  0.  This  illustrates  that  assumptions  A3 
and  A4  are  necessary  for  the  techniques  in  this  paper  to  be  directly  applicable.  An  alterna¬ 
tive  is  to  develop  a  theory  which  maintains  the  same  notion  of  admissibility,  but  relaxes  A3 
£md  A4.  The  easiest  way  to  do  this  would  be  to  pursue  the  suboptimal  case  introducing  e 
perturbations  so  that  A3  and  A4  are  satisfied. 

5.2.2  Relaxing  A1 

if  assumption  Al  is  violated,  then  it  is  obvious  that  no  admissible  controllers  exist.  Suppose 
A1  is  relaxed  to  allow  unstabilizable  and/or  undetectable  modes  on  the  ju  axis,  and  internal 
stability  is  also  relaxed  to  also  allow  closed-loop  ju  axis  poles,  but  A2-A4  is  still  satisfied.  It 
can  be  easily  shown  that  under  these  conditions  the  closed-loop  Woo  norm  cannot  be  made 
finite,  and  in  particular,  that  the  unstabilizable  and/or  undetectable  modes  on  the  ju  axis 
must  show  up  as  poles  in  the  closed-loop  system. 

5.2.3  Violating  Al  and  either  or  both  of  AS  and  A4 

Sensible  control  problems  can  be  posed  which  violate  Al  and  either  or  both  of  A3  and  A4. 
For  example,  cases  when  A  has  modes  at  s  =  0  which  are  unstabilizable  through  B2  and/or 
undetectable  through  C2  arise  when  an  integrator  is  included  in  a  weight  on  a  disturbance 
input  or  an  error  term.  In  these  cases,  either  A3  or  A4  are  also  violated,  or  the  closed-loop 
Woo  norm  cannot  be  made  finite.  In  many  applications  such  problems  can  be  reformulated 
so  that  the  integrator  occurs  inside  the  loop  (essentially  using  the  internal  model  principle), 
and  is  hence  detectable  and  stabilizable. 

An  alternative  approach  to  such  problems  which  could  potentially  avoid  the  problem 
reformulation  would  be  pursue  the  techniques  in  this  paper,  but  relax  internal  stability  to  the 
requirement  that  all  closed-loop  modes  be  in  the  closed  left  half  plane.  Clearly,  to  have  finite 
Woo  norm  these  closed-loop  modes  could  not  appear  as  poles  in  Ttw  The  formulae  given  in 
this  paper  will  often  yield  controllers  compatible  with  these  assumptions.  The  user  would 
then  have  to  decide  whether  closed-loop  poles  on  the  imaginary  axis  were  due  to  weights  and 
hence  acceptable  or  due  to  the  problem  being  poorly  posed  as  in  the  above  example. 

A  third  alternative  is  to  again  introduce  e  perturbations  so  that  Al,  A3  and  A4  are 
satisfied.  Roughly  speaking,  this  would  produce  sensible  answers  for  sensible  problems,  but 
the  behaviour  as  c  — »  0  could  be  problematic. 

5.2.4  Relaxing  A2 

In  the  cases  that  either  D\2  is  not  full  column  rank  or  i?2i  is  not  full  row  rank  then  improper 
controllers  can  give  bounded  Woo-norm  for  although  will  not  be  admissible  as  defined  in 
section  1.4.  Such  singular  filtering  and  control  problems  have  been  well-studied  in  W2  theory 
and  many  of  the  same  techniques  go  over  to  the  Woo-case  (e.g.  Willems(1981),  Willems  et 
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oi.(1986)  and  Hautus  and  Silvennan(1983)).  In  particular  the  structure  algorithm  of  Silver- 
man  (1969)  could  be  used  to  make  the  terms  Du  and  Z?2i  full  rank  by  the  introduction  of 
suitable  differentiators  in  the  controller. 

5.2. 5  Behaviour  at  s  =  oo 

It  has  been  assumed  in  Theorem  5.1  that 

and  a  necessary  condition  for  a  solution  is  that  this  holds  with  >.  H  equality  holds  then  one 
or  both  of  the  Hamiltonian  matrices  cannot  be  defined.  This  corresponds  to  the  case 

ini  a{J^i{G{co),K{(Xi)))=  1 

if(oo) 

where  iif(oo)  is  just  considered  to  be  an  arbitrary  matrix.  If 
inf  K{juj)))  <  1,  for  some  w  =  Uo 

K(ju) 

then  a  bilinear  tr2msformation  from  the  right  half  plane  to  the  right  half  plane  that  moves  the 
point  ju>o  to  oo  win  enable  the  Hamiltonians  to  be  defined.  One  of  them  will  however  have  an 
eigen  value  at  the  point  on  the  imaginary  axis  to  which  the  point  at  oo  has  been  transformed. 
A  more  intricate  situation  arises  when 

■K’O'u?)))  =1,  V  w. 

K(jw) 

Here  the  corresponding  J-factorization  problem  (see  Green  et  of. (1988))  or  spectral  factoriza¬ 
tion  problem  is  rsink  deficient  for  all  u;.  The  theory  of  spectral  factorization  for  such  cases  can 
be  derived  via  the  solutions  to  a  Linear  Matrix  Inequality  (Willems(1971)),  or  via  the  stable 
deflating  subspace  of  the  zero  pencil  (see  Van  Dooren  (1981)  and  Clements  and  Glover(1989)). 
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Introduction 


The  space  station  is  acted  upon  by  internal  torques  and  external  torques.  The  internal  torques  are 
due  to  moving  payloads,  astronauts,  control  moment  gyros  (CMGs),  etc.  The  external  torques  are  due 
to  aero,  gravity  gradient,  earth’s  magnetic  field,  reaction  jets,  etc. 

If  there  were  no  external  torques,  the  total  torque  on  the  system  could  be  kept  zero  by  command¬ 
ing  the  CMG  torque  to  offset  the  other  internal  torques.  The  resulting  CMG  momentum  would  remain 
finite  since  the  momentum  of  the  other  internal  torque  producing  elements  are  finite. 

To  keep  the  system  momentum  bounded,  the  average  external  torques  must  also  be  zero.  This 
could  be  done  by  using  the  reaction  Jets,  however  they  use  expendable  fuel.  In  order  to  minimize  the 
use  of  expendable  fuel,  some  other  source  of  external  torque  must  be  used  to  counter  the  aero  torque. 
The  large  moments-of-inertia  (MOI)  of  the  space  station  give  rise  to  attitude  dependent  gravity  gradient 
torques  large  enough  to  offset  the  aero  torque. 

The  gravity  gradient  torques  depend  on  the  MOI  of  the  vehicle,  so  variations  of  the  MOI  can  have  a 
significant  impact  on  the  closed  loop  system  stability  and  performance. 

Several  earlier  studies  of  the  space  station  attitude  dynamics  have  looked  at  ways  of  analyzing  the 
mass  properties  variations  which  appear  in  the  moment-of-inertia  matrix  [WBWGLS],  [BWGS],  [WW], 
[RS],  [BP],  Since  the  variations  are  quite  large,  it  is  important  to  remove  as  much  conservatism  as 
possible  by  taking  advantage  of  the  known  structure  of  the  moment-of-inertia  matrix.  The  structure  of 
a  generic  moment-of-inertia  matrix  is  completely  described  by  the  fact  that  it  is  symmetric  positive 
definite  and  the  sum  of  any  two  of  its  eigenvalues  is  equal  to  or  greater  than  the  remaining  eigenvalue. 
More  specific  structure  is  implied  if  the  mass-properties  variations  arise  from  some  specific  scenario 
such  as  a  single  payload  mass  moving  in  some  presctioed  bounded  region. 

In  this  paper,  three  cases  are  examined  in  detail.  The  first  case  exactly  represents  the  perturbed 
system  dynamics  for  all  moment-of-inertia  variations  due  to  a  payload  of  fixed  mass  moving  in  a  given 
rectangular  region.  The  second  case  exactly  represents  the  perturbed  system  dynamics  for  all  possible 
moment-of-inertia  matrices.  The  third  case  exactly  represents  the  perturbed  system  dynamics  for  all 
diagonal  moment-of-inertia  matrices.  Case  1  can  be  combined  with  either  case  2  or  case  3. 
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Expressing  the  Mass-Properties  Perturbations  as  Linear  Fractional  Transformations 

We  have  been  using  |x  (stnictured  singular  values)  to  analyze  the  robustness  of  the  space  station 
controllers.  In  order  to  use  the  p-synthesis  or  p-analysis  techniques  on  a  perturbed  linear  system,  the 
perturbations  must  appear  as  rational  functions  of  some  undetermined  parameters  5.  Then  linear  frac¬ 
tional  transfonnations  can  be  used  to  put  the  perturbed  system  into  the  standard  A,  P,  K  format  used 
in  the  \i  computations  [ZD],  This  procedure  is  related  to  the  factorization  of  matrix  polynomials  in 
many  5  variables  and  can  lead  to  very  large  dimensions  of  the  A  matrix  (especially  if  the  factorization 
is  not  minimal). 

The  linearized  angular  dynamics  of  a  rigid  body  are  rational  functions  of  the  moment-of-inertia 
matrix  J,  so  if  the  perturbed  J  can  be  written  as  a  rational  function,  then  the  perturbed  dynamics  can  be 
put  into  the  form  required  for  the  |X  computations. 

The  variation  in  the  moment-of-inertia  matrix  J  can  come  from  several  different  sources.  The  pri¬ 
mary  variations  in  J  come  from  movable  payload  masses  and  from  variations  in  the  mass  properties  of 
the  space  station  itself. 

Three  types  of  variations  in  J  are  considered.  Case  1  considers  variations  due  to  a  movable  pay- 
load.  This  type  of  variation  is  quadratic  in  the  payload  location  and  linear  in  the  payload  mass,  so  they 
come  naturally  in  the  form  required  for  p.  computations.  A  detailed  description  is  given  for  how  to 
represent  this  perturbed  system  in  the  A,  P,  K  format.  Care  was  taken  to  find  the  smalles'  possible 
dimension  for  the  required  A  matrix,  however  the  resulting  A  matrix  was  still  30  x  30.  This  perturba¬ 
tion  structure  was  used  to  compute  p  for  various  space  station  attitude  controllers  and  produced  reason¬ 
able  results. 

Case  2  considers  generic  variation  in  the  J  matrix.  Let  J  be  Represented  in  a  factored  form: 
J  =  U  r  U^,  where  U  is  a  special  (det  =  1)  orthogonal  matrix  and  Z  is  a  positive  diagonal  matrix.  The 
physical  origin  of  the  J  matrix  ensures  that  the  entries  of  the  L  matrix  are  linear  in  the  total  mass  and 
quadratic  in  the  geometric  distribution  of  the  mass.  In  order  to  express  the  entire  J  matrix  as  a  rational 
function,  we  must  also  be  able  to  express  the  U  matrix  as  a  rational  function  of  its  three  independent 
parameters.  This  is  made  possible  by  the  Cayley  Transform.  Details  of  representing  generic  variations 
in  J  as  rational  functions  are  given.  The  smallest  size  A  matrix  for  this  generic  case  is  approximately 
160  X  160  which  we  considered  to  be  too  large  for  practical  use.  However  if  the  J  matrix  is  assumed 
to  be  diagonal,  then  the  problem  simplifies  enough  to  be  of  practical  use. 

Case  3  considers  all  possible  diagonal  perturbations  to  the  J  matrix. 


3 


NOTATION 

MASS  PROPERTIES 

X 

r  =  y  The  vector  from  the  space  station  c.g.  to  the  payload  c.g. 


0  -z 

y 

t  =  z  0 

-X 

.-y  X 

0. 

mi 

The  mass  of  the  space  station  (without  payload) 

m2 

The  mass  of  the  payload 

m 

The  reduced  mass  of  the  space  station  and  payload 

Jj  The  moment-of-ineitia  matrix  of  the  space  station  (without  payload) 

Jj  The  moment-of-inertia  matrix  of  the  payload 

Ji2  =  -m?^  The  increase  in  MOI  due  to  separation  of  the  space-station  and  payload  c.g.’s 
J  The  system  moment-of-inertia  matrix  with  respect  to  the  system  c.g. 

^xy  ^xz 

J,y  Jyy  Jyj  Thc  componcnts  of  the  system  moment-of-inertia  matrix 

■^xz  ^yz  ^zz 

STATES  AND  INPUTS 

All  veaors  without  superscripts  are  in  body  axes.  Vectors  with  an  LV  superscript  are  in  the  Local 
Vertical  Local  Horizontal  reference  frame  which  is  centered  in  the  spacecraft  and  rotates  at  orbit  rate  (x 
axis  along  the  flight  path,  z  axis  towards  the  earth,  y  axis  perpendicular  to  the  orbit  plane), 
m  angular  rate  vector 

C  =  [^  Cy  rotation  matrix  from  the  LVLH  reference  frame  to  the  body  axes  reference  frame. 

H  CMG  momentum  vector 

Imio.  Zcmg  Aero  and  CMG  torques 

LINEARIZED  STATES  AND  INPUTS 

At  equilibrium,  the  principle  body  axes  will  not  be  aligned  with  LVLH,  so  define  a  new  set  of  body 
axes  which  are  aligned  with  LVLH  at  equilibrium.  The  MOI  matrix,  J,  will  not  be  diagonal  in  these 
new  body  axes.  The  set  of  three  angles,  0  are  the  small  deviations  from  equilibrium. 

0 

0  small  angle  deviations  between  the  LVLH  and  body  reference  frames  (not  a  vector) 

h  linearized  momentum 

Ixero*  IcMG  linearized  aero  and  CMG  torques 


MISCELLA  jEOUS 


I 


n 


LV 


0 

-Oq 

0 


ei  = 


1 

0 

0 


62  = 


0 

1 

0 


63  = 


0 

0 

1 


An  n  X  n  i(l6ntity  matrix 

angular  rat6  of  th6  LVLH  r6f6r6nc6  fram6  (orbital  rat6) 
a  unit  v6Ctor  along  th6  flight  path  in  LVLH  coordinat6S 

a  unit  v6Ctor  normal  to  th6  orbit  plane  in  LVLH  coordinates 

a  unit  vector  pointing  towards  earth  in  the  LVLH  reference  frame 
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Nonlinear  Equations  of  Motion 


The  nonlinear  equations  of  motion  for  the  space  station  and  a  moving  payload  have  been  examined 
in  [LSS]  and  [WHS],  These  references  examined  the  results  of  various  payload  motions  using  noii- 
linear  simulations.  In  this  discussion,  we  will  assume  that  the  speeds  and  accelerations  of  the  payload 
are  small,  however  the  separation  of  the  space  station  c.g.  and  the  payload  c.g.  can  be  quite  large  so 
the  resulting  gravity  gradient  torques  will  be  quite  large.  In  later  sections,  we  will  show  how  to 
represent  the  linearized  dynamics  as  a  linear  fractional  transformation  in  the  (x,y,z)  position  of  the  pay- 
load.  This  will  allow  us  to  use  structured  singular  values  to  determine  how  robust  a  control  system  is 
to  the  MOI  variations. 


diagram  of  earth  and  2-body  SSF/payload 


If  the  equations  of  motion  are  written  with  respect  to  the  v.omposite  system  c.g.,  then  the  transla¬ 
tional  and  rotational  dynamics  decouple.  The  equations  are  further  simplified  by  assuming  that  the  two 
bodies  have  the  same  angular  velocity. 

The  torque  due  to  gravity  gradient  is  given  by: 

Igg  =  3<Oo^  6*  J  e^  (1) 

The  torque  due  to  payload  velocity  and  acceleration  is  given  by: 
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(2) 


T^y  =  mT  2a  r  - r 

We  will  assume  that  payload  motion  with  respect  to  the  main  body  is  usually  quite  slow,  so  that  Tj.” 
can  be  represented  as  a  bounded  external  disturbance.  However  the  MOI  matri\,  J,  depends  on  the 
payload  position  r  and  the  resulting  changes  in  the  gravity  gradient  torques  can  be  quite  large. 

We  will  assume  that  the  aero  torque  can  be  represented  as  a  bounded  external  disturbance. 

The  attitude  dynamics  are  give'  by: 

J  m  +  &  J  ©  =  3©o^  £z  “  Z 

where 

Z  =  ZcMG  “  Z^.V  ~  Zaero  (4) 

The  attitude  kinematics  are  given  by: 

C  +  (a  -  %)  C  =  0  where  C  =  [e^j,  Cy,  ej  (5) 

The  CMG  dynamics  are  given  by: 

H  +  a  H  =  ZcMG  (6) 
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Linearized  Equations  of  Motion 


Note  that  at  equilibrium,  the  principle  body  axes  can  have  arbitrary  orientation  with  respect  to  the 
LVLH  reference  frame,  so  we  will  define  a  new  set  of  body  axes  which  are  aligned  with  LVLH  at 
equilibrium.  In  tliis  new  set  of  axes,  C  =  1  at  equilibrium,  but  J  is  not  diagonal. 

Since  C=I  at  equilibrium,  %  =  -coq^,  and  =  63. 

The  equilibrium  equations  are: 

0  =  — ^  Je2  +  3(Oo^  63  “  Xo  (7) 

solving  for  the  equilibrium  torque,  we  get: 


We  will  now  determine  the  linearized  equations  for  small  variations  about  this  equilibrium. 
For  first  order  variations  from  equilibrium, 

C  =  I  +  0  where  8  is  skew  symmetric 

C=  ^0 
dt“ 


(9) 


Substituting  this  into  the  nonlinear  attitude  kinematics,  and  keeping  first  order  terms  in  0  and  0  gives: 

®  =  0  +  0  (10) 
Differentiating  this  and  keeping  first  order  terms  in  0  and  0  gives: 

m  =  0  +  ^^''0  (11) 


We  now  have  expressions  for  the  first  order  variations  of  the  states  and  their  derivatives  in  terms  of  the 
angles  and  rates  0,  0  between  the  LVLH  reference  frame  and  the  body  reference  frame.  The  torque 
inputs  must  also  be  written  as  first  order  variations  from  the  equilibrium  torque,  however,  the  results 
depend  on  whether  we  use  the  body  reference  frame  or  the  LVLH  reference  frame  to  express  these 
torques.  Note  that  Tq  =  Tq^'^. 

If  we  use  the  body  reference  frame  for  the  torques  (as  in  [WBWGLS]),  we  get 

Z  =  Xo  +  ^3^G  (12) 

If  we  use  the  LVLH  reference  frame  for  the  torques,  we  get: 

T  =  C  T^v  =  (I  +  5)(ToLV  ^  ^^LV)  ,  To  +  Tcmg"'''  -  Xo  §  (13) 

Since  the  -Tq  term  multiplies  the  first  order  term  0,  it  must  be  included  in  the  A  matrix  of  the  linear¬ 
ized  system  when  using  the  LVLH  torque  inputs. 
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The  equilibrium  value  of  the  torque  is  given  by: 


lo  =  <Bo^(363  J  63  -  §2  J  62)  =  ©0^  3J. 


^  2 


Note  that 


=  €■*■  a  C 


Plugging  the  first  order  terms  for  the  states  and  the  inputs  into  the  nonlinear  attitude  dynamic  equations 
and  keeping  only  the  first  order  terms  in  0,  0,  0  and  gives: 

J  (  0  +  =  -  0  +  §]  X  [  J  [0  +  6  + 

f  irr  Ilf  1 

3(0o^  £3  “  |3§  X  J  ^  - 130  -  To  + 1 

L  JLL  *^JL  J 

=  X  [  J  I  +  3(0n^  e^  X  r  J  e^  I  -  To  constant  terms  n 


=  -©0^'^  X  J  +  3(Oo^  63  X  J  63  j  -  To  const 

+  -mo^'^  X  J  !§  +  ao^''  0  +  3ao^  ^  x  J  |^  -  ^0 

+  J  Wo^''  X  0  +  0  -  3(Oo^  "I 


linear  terms 


linear  terms 


+  higher  order  terms 


Keeping  only  the  linear  terms  gives: 


IJ  0  0  0 


I  -I2 


r  1 0  J  0  0  0  -53  r^j  0 

'<2-“»529=[la  %-I  3-JoO  Ji,  0  I  -S,  Ji 


§2  ®0  § 


1 0  0  0  Je3  0  63 
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J  0  =  [Le(J).  Led)] 

=  [Le(J).  Led)] 

=  [Le^^'d).  LeLV(j)] 


©0  0 
©0^0 


~  ICMG 


“00  _r 
©0^  0  I- 


ICMG 


LV 


-ioe] 


©00 

©0^0 


“ICMG 


LV 


(16c) 


where  and  L(J)  are  the  following  linear  functions  of  J. 


L^^'d)  =  [Le^^d).  Le^'^d)]  = 


[Jh-  +  I2J  .  3^[J  63  -  (J'ej)]  -yiJh-  (J’ez)] 


0 

-21 


yz 


2Jyz  Qxn  lyy  Izz) 
0  2J 


(-J«  +  Jyy-Jzz)  -2J 


xy 


'xy 

0 


^(Jyy  -  41 


'yy 

3J 

-J, 


xy 


xy 


-3d„  -  U  -3J 


'yz 


yz 

Ixx  ~  •^yy 


(17) 


L(J)  =  [Led).  Led)] 


where 


Led)  =  L0^'^(J)  skew  symmetric. 

IT 

■=|Le"''d) 

©o‘ 


Led)  =  Le^'^d)  -  -%  =  L^'^d) 


(18) 


(19) 


In  the  nonlinear  equations  of  motion,  the  gravity  gradient  torques  were  not  a  function  of  the  rotation 
angle  about  the  e^  vector  and  the  gravity  gradient  torques  had  zero  component  in  the  z  direction.  In 
the  linearized  equations  of  motion  we  get  one  or  the  other  of  these  attributes,  depending  on  whether  we 
represent  the  torques  in  the  LVLH  or  body  axes  reference  frame. 

For  inputs,  the  gravity  gradient  torque  has  zero  component  about  the  z  axis  which  is 

aligned  with  the  gravitational  force.  However,  as  the  0^  angle  varies,  the  gravity  gradient  torque 
rotates  around  in  the  x,y  plane. 

For  Xcmg  inputs  the  gravity  gradient  torques  are  not  a  function  of  0^,  but  they  have  components 
in  all  three  axis. 

Linearization  of  the  CMC  dynamics  is  simpler  in  the  LVLH  reference  frame: 

=  +  (20) 
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diagram  of  J  and  L(J)  and  integrators 


11 


Moments-of-Inertia  for  the  Composite  Vehicle 


The  system  eg  is  given  by: 


leg  ~ 


_  jjgl  +  ^2  leg?. 


mi  +  m2 


(21) 


and  the  MOI  of  the  system  about  the  system  eg  is  given  by: 

J  =  -  J 


system 


=  -/(&-  l,,)^()m{R)  -  J  (R  -  -  !«i)' 

body]  body2  ^2 


(22) 


Let 


Ji=-  f  (R-T^gi)^dm® 

bodyl 


h  =  -  j  (R-l:g2)^cim(R) 

body2 


(23) 


and 


J12  =  - 


mi  m2 


- —^<le»2  -  Icpl) 

mi  +  m2 


(24) 


then 


J  —  Ji  +  J2  +  J12 


(25) 


We  ean  simplify  the  expression  for  J12  by  setting 


m  = 


mj  m2 
mi  +  m2 


aiKl 


*"  *’cg2  ^cgl 


givmg 


J12  =  -m  = 


m 


y^+z^ 


xy  -xz 
-xy  x^+z^  -yz 
-xz  -yz  x^+y^ 


(26) 


(27) 
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Perturbed  System  Dynamics  due  to  Mass-Properties  Variations 

The  only  dependence  on  J  in  the  linearized  dynamics  comes  from  the  equation  for  6. 

We  will  now  show  how  variations  in  Jj,  J2,  and  J12  can  be  combined. 

Since  L(J)  and  are  linear  in  J, 

L(J)  =  L(J,)  +  UW  +  L(Ji2) 
and 


so  the  linearized  dynamics  can  be  written  as 


(Ji  +  J2  +  J12)  e  =  L(Ji)  +  L(J2)  +  L(J,2) 


a>o9 

0)0^  e 


ICMG 


or 


(Ji  +  J2  +  J12)  e  =  +  L^''(Ji2) 


©0  e 

©0^  e 


“ICMG 


LV 


(28) 


(29) 


diagram  of  Ji  J2  J12  UJi)  L(J2)  L(Ji2)  and  integrators 


We  will  assume  that  m,,  m2,  J2.  and  ©0  are  all  known  constants.  We  will  look  at  the  following  three 
cases  for  J}  and  Ji2- 
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CASE  1:  J]  is  a  known  constant,  but  J12  is  perturbed  (variable  payload  position). 

CASE  2:  J12  is  a  known  constant,  but  Jj  has  arbitrary  perturbations  (variable  core  body  MOl). 

CASE  3;  J 12  is  a  known  constant,  but  Jj  has  diagonal  perturbations  (variable  core  body  MOl). 

Note  that  case  1  can  be  combined  with  case  2  or  case  3. 

Case  1:  Moment-of-Inertia  Variations  due  to  Payload  Position 

Putting  the  MOl  and  Torque  Matrices  Directly  into  Linear  Fractional  Form 

J12  has  only  three  independent  parameters,  (x,y,z),  so  we  must  write  J12  and  L(Ji2)  as  linear  frac¬ 
tional  transformations  with  respect  to 

Kinx 

A  =  5y  (30) 

5,1^ 

where  Inx,  I„y,  and  are  identities  of  the  appropriate  dimensions. 

The  first  step  is  to  write  J12  and  L(Ji2)  as  explicit  functions  of  (x,y,z).  Equation  27  already  gives  J12 
as  an  explicit  quadratic  function  of  (x,y,z). 

Using  Equation  27  in  equation  17,  we  can  get  L^\Ji2)  as  an  explicit  quadratic  function  of  x,  y,  z. 

0  -2yz  2y^  4(y^  -  z^)  -4xy  4xz 

11^)  =  m  2yz  0  -2xy  ,  m  -3xy  3(x2  -  z^)  3yz  (31) 

-2y^  2xy  0  xz  -yz  y^  -  x^ 


Polynomial  (or  even  rational)  matrix  functions  in  several  variables  can  be  put  directly  into  the 
linear  fractional  format,  however  the  dimension  of  the  required  A  can  be  very  large  [ZD]. 

In  order  to  write  as  a  linear  fractional  transformation  with  respect  to  A,  first  write  r  as  an 
average  value  plus  differences  times  Sj. 

Assume  that  the  payload  is  restricted  to  move  in  a  rectangular  region  defined  by  the  following  set: 
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• 

XdifSx 

15,1  <  1.0 

re  ■ 

Live 

ydif  Sy 

with  I5yl  <  1.0 

Zdif  5z 

I5J  <  1.0 

• 

(32) 


The  3x6  matrix  polynomial  in  equation  31  can  be  written  using  3x6  coefficient  matrices 
Qij.  Qi.  Qo  which  are  ftmctions  of  r^*  and  and  z^. 

i=l 


SQij5d+Qi' 


5i  +  Qo 


(33) 


This  matrix  polynomial  can  be  written  as  a  linear  fractional  transformation  as  follows: 

[aq  BqI 


=  LFT,, 


Co  Df 


.Aq 


—  Dq  +  CqAq(I  —  AqAq)  *Bq 


where  Aq  is  chosen  to  be  rulpotent  of  order  1,  so  that  (I  -  Aq/^q)  *  =  (I  +  AqAq) 


(34) 


(4x3 

Qll 

(4x3 

Qi2 

(4x3 

Qi3 

(^6x3 

(^6x6 

(^6x3 

(^6x6 

(^6x3 

06x6 

(4x3 

(4x6 

(4x3 

Q22 

(4x3 

Q23 

(^6x3 

(^6x6 

(^6x3 

(^6x6 

(^6x3 

06x6 

(^3x3 

(4x6 

(4x3 

(4x6 

(4x3 

Q33 

(^6x3 

(^6x6 

(^6x3 

(^6x6 

(^6x3 

06x6 

Bq  = 


^3x6 

l6  ' 

03x6 

4x6 

O3X6 

4x6 


Cq=[i 


4x3  Ql  4x3  Q2  4x3  Q3 


-  Qo 


(35) 


(36) 


A  =  diag(5,l9.  5yl9.  5,19)  (37) 

This  representation  of  requires  27  5s.  In  addition,  it  takes  12  5s  to  represent  J  (or  r*)  as  an 
LFT  (see  the  section  on  expressing  T  as  an  LFT).  The  total  number  of  5s  required  is  39. 


In  order  to  reduce  the  dimension  of  the  A,  it  is  sometimes  useful  to  first  factor  the  matrix  polynomials. 
Polynomial  functions  in  several  variables  cannot  always  be  factored  into  linear  terms,  however  in  our 
case  it  is  possible. 
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Factoring  the  MOI  and  Torque  Matrices 


Ji2  is  already  in  faaored  form  (see  equation  27).  The  3x6  matrix  L(-[VmTj^)  can  also  be  factored 
into  left  and  right  factors,  each  of  which  is  linear  in  Vinr.  The  factors  are  not  unique.  The  freedom  in 
the  factorization  is  parameterized  by  an  arbitrary  vector  X  whose  values  are  then  chosen  to  reduce  the 
number  of  5s  needed  in  the  LFT  of  the  factor  (see  an)endix  for  details). 

L(-[VmT]2)  =  Frjyl^^ 

L^''(-[>/mT]2)  =  FL^'^(>/mi)  FR^''(V5ii) 


where 


FL('^r)  =  T 


^  _  fly  0  0  0  -3z  -yl  r  -i 

FR(''rai)  =  "Vm  0  2y  0  4z  0  -x  +'^  rXA^,  rXe^ 
LO  0  2y  4y  -3x  Oj  •-  -I 

r2y  0  0  0  — 4z  — 4yl  r 

Fl^''(VS[)  =  Vm  0  2y  0  3z  0  3x 

[O  0  2y  y  X  0  J  L"  -  -  - 

_  I  V/  I —  \  ^  ^3x3 

0,.; 


choose  Xq  =  0  =  Xe 


choose  Xq  = 


—  0 


Equation  16  now  becomes: 


[j,  +  Jj  -  (VS^j  0  =  [l(J,)  +  L(J2)  +  Fl(VEi)  FR(VSr)]  -  Tcmg 

=  [l'-'-CJ,)  +  K?j  - 


The  only  dependence  of  the  dynamics  on  m  and  r  is  through  Vin  T  and  F(Vm  r),  both  of  which  are 
linear  in  >^nr. 
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big  diagram  with  Jj  J2 
L(Ji)  L(J2)  Fsys  integrators 

for  body  and  LVLH  torque  inputs 


with  DELTA  sizes  labeled 


Putting  the  Factors  into  Linear  Fractional  Form 
Expressing  as  a  Linear  Fractional  TVansformation 


So 

VmT  =  qA.,  Bj  +  D- 

where 


0  0  y^if  0  zjif  0 

^  Xjif  0  0  0  0  Zjif  Dl  =  Vin  ^ 

0  JW  0  ydif  0  0 

and 

Af  =  diag[5,l2,  Sylj, 


0  0-1 
0  1  0 
0  0  1 
-I  0  0 

0-10 
1  0  0 


(44) 


(45) 


(46) 


Let  be  the  following  "system"  matrix 


(47) 


Then  equation  44  becomes 

V5iT  =  LFTuj^P^,  A.^  (48) 

where  LFTu  is  a  Linear  Fractional  Transformation 

Since  J,2  =  the  above  LFT  for  T  is  used  twice  in  the  expression  for  for  J,2  on  the  left  hand 

side  of  equation  29. 
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Expressing  as  a  Linear  Fractional  Transformation 


Fr(V5i  r)  =  FR(Vm  r,ve)  + 


Vm  Xdif 


^  Ydif 


^  Zjif 


0  0 
1  0 
0  1 


5x  0 

0  5, 


0  0  1 


1  0 
0  1 
0  0 


0  0  0  0  0  -1 
0  0  0  0  -3  0 


2  0  0  0  0  -1 
0  2  0  0  0  0 
0  0  2  4  0  0 


8y  0  0 

I 

0  5y  0 

1 

- 

0  0  6y 

1 

5z 

0 


0  0  0  0  -3 
0  0  0  4  0 


(49) 


so 

Fr(aG  r)  =  Cp|j  Ap^  Bpij  +  Dp^  (50) 

where 


Br 


0  0  0  0  0  -1 
0  0  0  0  -3  0 
2  0  0  0  0  -1 
0  2  0  0  0  0 
0  0  2  4  0  0 
0  0  0  0  -3  0 
0  0  0  4  0  0 


Cp-  =  Vm 


0  0 
Xdif  0 

0  xjjf 


ydif 

0 

0 


0  0  Zjjf  0 

Ydif  0  0  Zdif 

0  ydif  0  0 


DF^  =  FR(VEr^ve)(51) 


and 


Ap^  =  diag[5,l2,  6yl3.6^2] 


(52) 


Let  Pp|^  be  the  following  "system"  matrix 


Pf,= 


0  Bp^ 
^Fr  ^Fr 


(53) 


Then  equation  50  becomes 

Fr(VS  r)  =  LFT„[Pf^.  Ap^ 


(54) 


where  LFTu  is  a  Linear  Fractional  Transformation 
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Expressing  as  a  Linear  Fractional  Transformation 


FL‘'''(VSr)  =  FL''''(>Sr,ve) 


+  \m  xjif 


+  Vin  ydif 


+  ^Zdif 


0  0 
1  0 
0  I 

’l  0 
0  0 
0  1, 

1  0 
0  1 
0  0 


5x  0 

0  5, 

5y  0 

0  Sy 


5^  0 
0  5, 


[-2  0  0  0  0  3 

LO  00010 

[2  0  0  0  0  -4' 

[0  0  2  1  0  0 


0  0  0  0 
0  0-23 


-4 

0 


(55) 


so 

Fl^'^(’S  r)  =  Cp^Lv  Ap^Lv  Bp^Lv  +  Dfl*-''  (56) 

where 


-2 

0 

2 

0 

0 

.0 


0  0  0 
0  0  0 
0  0  0 
0  2  1 
0  0  0 
0-2  3 


0  3 

1  0 
0  -A 
0  0 
-4  0 
0  0 


0 

Xdif 

0 


0 

0 

Xdif 


ydif  0  zjif  0 

0  0  0  Z(jjf 

0  ydif  0  0 


(57) 


Ap^Lv  =  diag[5,l2,  5yl2,  6^2] 


(58) 


Let  Pp^Lv  be  the  following  "system"  matrix 


0 

Cp^LV  Dp^LV 


(59) 


Then  equation  56  becomes 


FL^''(VSr)  =  LFT„ 


(60) 


where  LFT„  is  a  Linear  Fractional  Transformation 
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Case  2:  Generic  Mass-Properties  Variations  for  a  Single  Body 

A  generic  moment-of-inertia  matrix  can  be  written  as 

myy  +  mzz  0  0 

J  =  U  0  mxx  +  mzz  0  (61) 

0  0  mxx  -I-  myy. 

where  U  is  an  orthogonal  matrix  (  three  free  parameters  in  U  )  Note  that  this  parameterization  of  J 
automatically  satisfies  the  constraint  that  the  sum  of  any  two  of  its  eigenvalues  is  greater  or  equal  to 
the  third  eigenvalue. 


The  Cayley  Transform  provides  a  way  of  parameterizing  U  as  a  rational  function  of  three  parameters. 
Since  rational  functions  can  be  written  as  linear  fractional  transformations,  they  can  be  handled  by  the 
structured  singular  value  algorithms. 


Let  ^  e  and  let 


u  =  as  +  Dds  - 


(62) 


Then  U  is  an  orthogonal  matrix,  (U^  =  ir‘). 

Proof: 

(I3  +  1)  [uT"  -  U-‘]  (I3  1) 

=  (I3  +  D  [ih  +  (I3  -  D  -  (h  -  Dds  +  Ir]  (I3  + 1)  ^^3) 

=  (I3  -  ^)  -  (I3  -  ^)  =  0 

and 

det(l3  +  ^  =  l-K  1^1  ^  0  (64) 


so 

ir‘  -  =  0  (65) 

The  LFTs  for  (I3  +  ^)  and  Os  -  ^)'’'  each  take  6  8s,  so  U  and  each  take  12  deltas.  Z  also 
takes  6  6s,  so  J  takes  a  total  of  30  5s.  Since  the  expression  for  the  dynamics  in  equation  16b  has  J  in 
it  5  different  places,  it  could  take  as  many  as  as  150  &  for  the  whole  system  when  generic  moment- 
of-inertia  matrices  are  used. 
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Case  3:  Diagonal  Mass-Properties  Perturbations  for  a  Single  Body 

In  [BP],  MOI  variations  which  leave  J  diagonal  arc  put  into  LFT  form.  Some  of  those  results  arc 
reproduced  here. 

If  U  =  I,  so  J  is  always  diagonal,  then  the  moment-of-inertia  matrix  can  be  written  as 


J  = 


myy  +  mzz  0  0 

0  mxx  +  mzz  0 
0  0  mxx  +  myy 


(66) 


Note  that 


Jyy  +  1x2  JjUt 

mxx  =  — - - - 

2 


myy  = 


fxx  fz2  ■^yy 


mzz  = 


■^xx  Jyy  Jzz 


(67) 


For  this  J. 


lLV(J)  =  L(J)  = 


0 

0 


0  2myy 
0  0 


-2myy  0  0  J 


4(rayy  -  mzz)  0  0 

0  3(mxx  -  mzz)  0 

0  0  myy  -  mxxj 


(68) 


Since  J  and  L(J)  are  both  linear  in  the  three  independent  parameters  mxx,  myy,  and  mzz,  they  can  be 
put  directly  into  linear  fractional  form. 


Assume  that  each  diag  element  of  J  lies  between  some  maximum  and  some  minimum. 
Let 


Jave  “ 


Jmax  Jmin 


Jdif  “ 


Jmax  Jmin 


(69) 


Let 


mxxdif  = 

myydif  = 

mzzjif  = 


Jdifyy  +  Jdif„  -  Jdif„ 

2 

2 

Jdif„  +  Jdify^  -  Jdif„ 

2 


(70) 
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where 


0  d 
1  0 
P  1. 

1  C 
0  C 
.0  1. 

1  d 
0  1 
P  Q 


^nua 

0  ■ 

0 

mxxjif 

0 

0 

5mM 

0 

0 

mxxaif 

^yy 

0  ■ 

[myydif  0 

0 

0 

^myy 

i 

0  0 

myydif 

+ 


+ 


5mzz 

0  ■ 

mzzdif 

0 

d 

0 

Smzz 

0 

mzZjjif 

0 

0 

Bj" 

■■ 

=  LFT„ 

.  Aj 
. 

(71) 


0 

mxxjif 

0 

0 

0 

mxx^ 

myydif 

0 

0 

0 

0 

1 

0 

1 

d 

0 

0 

myydif 

Ci  = 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

0 

mzzjjf 

0 

J 

0 

mzzjif 

0 

(72) 


Aj  diag(5njjjjjl2,  ^myy^2>  ^mzz^2) 


(73) 


L(J)  = 


0  0  2myy  4(myy  -  mzz)  0  0 

0  0  0  0  3(mxx  -  mzz)  0 

-2myy  0  0  0  0  myy  -  mxx 


=  L(J,ve)  + 


0  o' 

r 

1  0 

p  Ij 

“ 

1  b 

0  0 

P  1. 

“ 

1  b 

0  1 

P  0 

“'myy 

0 


0 

5n»x 

0  ‘ 
^myy_ 

0 

5mzz 


=  lft„ 


0  0  0  0  3mxXdjf  0 
0  0  0  0  0  -mxxjif 


0  0  2myydif  4myydif  0  0 

-2myydif  0  0  0  0  rnyy^jf 


y  L 


0  0  0  -4mzZdif  0  0 

J  0  0  0  -3mzzjif  0 


0  Bl(j) 

L^uj) 


.A, 


‘L(D 


(74) 


where 


0  0 

0  0 

0  0 

-2myydif  0 

0  0 

0  0 


0 

0 

Smxxdjf 

0 

0 

0 

0 

-mxxdif 

2myydif 

4myydif 

0 

0 

0 

0 

0 

myydif 

0 

-4mzz<iif 

0 

0 

0 

0 

-3mzzdif 

0 

Cud  - 


b  0  1  0  1  0 
1  0  0  0  0  1 
p  1  0  1  0  0 


(75) 


I7l(J)  -  L(Jave)  Al(j)  -  diag(5n^l2,  6n,yyl2,  5^12) 


(76) 


This  requires  6  6s  for  L(J)  and  6  6s  for  J. 
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Scaling  the  Inputs,  Outputs,  States,  and  Time 


Using  English  units,  the  space  station  has  moments  of  inertia  on  the  order  of  HJH  =  10*  and  the  orbi¬ 
tal  rate  is  coq  =  10~^.  These  numbers  show  up  raised  to  various  powers  in  the  state  space  equations 
and  give  them  large  condition  number.  When  solving  the  Ricatti  equations  for  the  or  H2  controll¬ 
ers,  numerical  problems  arise  due  to  the  poor  condition  number  of  the  state  space  equations.  One  way 
to  dramatically  improve  the  condition  number  is  to  make  the  system  dimensionless.  We  made  the 
equations  of  motions  dimensionless,  by  substituting; 


<flo— 4  1 


and 


J 


(77) 


(78) 


In  order  to  make  the  perturbation  equations  dimensionless  too,  we  must  also  substitute: 


This  is  equivalent  to  scaling  time,  scaling  the  state,  and  scaling  the  inputs. 


VARIABLE  SCALE  FACTOR 


time  l/(0o 

momentum 

angular  rate  O/g 

angles  1 

torques  Jzz®o^ 


(79) 
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Appendix  on  Factorization  of  3  x  3  Quadratic  Matrix  Polynomials 

Generic  3x3  homogeneous  quadratic  matrix  polynomials,  Q(r),  in  3  variables  rj,  T2,  may  not  be 
factorable,  however  we  are  concerned  with  a  special  class  of  such  matrices  which  satisfy  Q(r)  r  =  0. 
This  class  of  matrices  does  factor  into  left  and  right  linear  factor  matrices:  Q(r)  =  FL(r)  FR(r),  and  an 
explicit  formula  is  given  for  each  factor. 

THEOREM: 

Given: 

r  e  (3  variables) 

T  e  skew  matrix  formed  from  elements  of  r 

aiiv  9ii3>  922,.  g223>  333,.  9332.  923,.  9i3j  e  R^^  (24  coefficients) 

X  e  R^  (3  coefficients) 

F(r)  G  R^’^  Matrix  whose  elements  are  homogeneous  linear  polynomials  in  the  elements  of  r 

Q(r)  e  R^^  Matrix  whose  elements  are  homoget^ous  quadratic  polynomials  in  the  elements  of  r 

The  following  three  statements  are  equivalent: 

1)  Q(r)r  =  0 

2)  Q(r)  =  F(r)T 

3)  QQ)  has  24  free  coefficients  q,  and  is  of  the  form: 

Q(r)  =  0,  qiij.  91I3  922,.  Q.  922,  r2r2  +  j^933,.  933j.  0  r3r3  + 

+  923,.  -9223.  -933j  +  -9II3’  3132’  ~933,j  +  -^U^.  -922,-  “923,  “  9l32  i'ir2 

PROOF 
2)  ->  1) 

Q®  r  =  F(r)Tr  =  g 


3)  ->  2) 

Given  any  matrix  of  the  form  of  statement  3),  set  it  equal  to  a  general  3x3  homogeneous  linear 
matrix  times  T  and  then  equate  coefficients.  A  general  homogeneous  linear  matrix  polynomial  in  the 
entries  of  r  is  of  the  form: 
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F(r)  =  Fi  ri  +  Fj  rj  +  F3  Tj 
where  Fj  e  are  coefficient  matrices.  Setting 

I^Fi  ri  +  F2  rj  +  F3  r3  T=Q(r) 

and  equating  coefficients  gives: 

Fi  =  ^  -  §132>  “Sills’ 91I2  ^2  =  1:9223’ ^  +  923i’ “922,  =  “9332*  933,’ ^  +  9l32 

where  X  is  an  arbitrary  vector.  Note  that  F(r)  has  27  free  parameters,  q,^  so  aU  possible  homo¬ 
geneous  linear  matrices  can  be  generated  as  factors  of  the  quadratic  matrices  given  in  statement  3). 

1)  ->  3) 

A  general  3x3  matrix  homogeneous  quadratic  polynomial  in  three  variables  has  54  free  coefficients 
and  is  of  the  form: 

Q(£)  =  Qii  rirj  +  Q22  r2*‘2  +  Q33  r3r3  +  Q23  r2r3  +  Q13  rir3  +  Q12  rir2 

where  the  Qjj  are  3  x  3  coefficient  matrices. 

Let  the  three  columns  of  each  Qjj  be  denoted  by: 

Qij  =  9ij,’  9ij2’  9ij3 

then 

Q(^  r  =  gii^  rjrjri  +  3222  hhh  +  9333  ’■3r3’r3  + 

(922,  +  9122)  •'i*'2''2  +  (933,  +  9133)  rir3r3  +  (qii^  +  £12^)  ririr2  -1- 
(9332  +  9233)  ^2*’3r3  +  (91I3  +  913,)  +  (922,  +  9233)  ^^2^  + 

(923,  +  9i32  +  9123)  hhh 

Since  Q(r)  r  must  vanish  identically,  all  ten  vector  coefficients  must  be  zero,  leaving  only  8  free  vec¬ 
tor  quantities  to  parameterize  Q.  The  8  free  vector  quantities  are  those  in  the  expression  of  statement 
3). 
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The  constructive  formulas  for  parameterizing  Q(r)  and  F(i)  give  a  splitting  for  the  following 
homological-algebra  exact  sequences: 

Long  Exact  Sequence: 


The  kernel  of  each  map  is  exactly  the  range  of  the  previous  map. 
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ABSTRACT:  Stability  of  Dynamic  Inversion  Control  Laws  Applied  to 
Nonlinear  Aircraft  Pitch-Axis  Models 
by  Blaise  Morton  and  Dale  Enns,  Honeywell  SRC,  Minneapolis 


Introduction 


Dynamic  inversion  is  a  nonlinear  control  technique  that  has  been  ^lied  by  Honeywell  to  a  variety  of  realistic 
aerospace  vehicle  models  with  reasonably  good  results.  The  list  of  study  applications  includes  models  of  the  F- 
14  aircraft,  the  HARV  F-18  aircraft,  a  McDonnell  Douglas  model  of  the  NASP  vehicle,  and  a  General  Dynamics 
model  of  a  next-generatitm  booster  vehicle.  The  main  advantage  of  dynamic  inversion  over  more  conventional 
linear  control  techniques  is  its  tq>plicability  to  the  fuU  nonlinear  vehicle  models. 

The  theory  of  dynamic  inversion  (and  nonlinear  control  in  general)  is  not  well  understood.  Most  of  what  we 
know  about  dynamic  inversion  theory  is  summarized  in  the  references  [El],  [E2],  and  [MEHH]. 

This  note  describes  a  global  stability  result  for  dynamic  invasion  applied  to  nonlinear  aircraft  pitch-axis  models. 
The  point  is  to  examine  the  technique  from  a  mathematical  point  of  view  and  try  to  understand  why  and  how  it 
works. 


Section  1:  Equations  of  Motion 

We  concentrate  on  aircraft  pitch-axis  models  similar  to  those  used  in  current  aerospace  vehicle  design.  The 
body-axis  coordinate  system  is  used.  See  Figure  1.  There  are  four  states: 

U  =  component  of  velocity  in  the  aircraft  longitudinal  (x)  axis 
W  =  component  of  velocity  in  the  aircraft  vertical  (z)  axis 
Q  =  vehicle  pitch-rate 

e  =  vehicle  pitch  attitude  relative  to  local  horizontal. 

The  equations  of  motion  are: 


■  c,(a)  ' 

■  C,,(a)  ■ 

T_ 

m 

m 

U 

-WQ 

-g  sin(0) 

m 

c;(a) 

C.,6(a) 

d 

W 

UQ 

g  cos(0) 

0 

m 

1  -> 

m 

dt 

Q 

— 

0 

+ 

0 

+ 

0 

+  ipV  s 

cCkKo) 

+  ypV^S 

cCM.6(a) 

e 

Q 

0 

0 

Iy 

Iy 

0 

0 

The  variables  in  equation  (1.1)  have  the  following  meaning: 

g  =  gravitational  acceleration  (constant), 

T  =  thrust  (control  input), 

m  =  vehicle  mass  (constant), 

p  s  air  density  (^sumol  constant  hoe), 

V  =  speed  =  Vu^  + W®, 

w 

a  =  angle  of  attack  =  atan(-^), 
ly  =  vehicle  inertia  (constant). 
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S  =  elevator  angle  (control  input), 
c  =  mean  aerodynamic  chord  (constant). 

C,(a),  Ci(a).  CM(a)  =  aerodynamic  functions  for  5  =  0, 

Cx^(a),  Cx,5(a),  CM^(a)  =  aerodynamic  functions  due  to  nonzm  5. 

The  two  control  inputs  T  and  5  are  assumed  to  be  limited  to  values  within  a  fixed  int^al.  A  reasonable  set  of 
ranges  for  a  fighter  is  0  ^  T  ^  mg  and  -20  degrees  ^  5  <  20  degrees.  The  aerodynamic  functions  usually  depend 
on  Mach  as  well  as  a,  we  neglect  the  Mach  dqiendence  here.  Though  this  model  has  a  number  of  special 
features  introduced  to  simplify  analysis,  most  of  the  discussion  below  can  be  augmented  to  apply  to  models  of 
more  goieral  type. 


Section  2:  The  Equilibrium  Manifold 


The  model  (1.1)  has  the  following  form: 


X  =  f(x,u) 


(2.1) 


where  x  is  a  vector  in  R“  and  u  is  a  vector  in  R“.  Let  f»(x,u)  denote  the  Jacobian  matrix  of  the  transformation  f 
from  R“  X  R“  to  R“.  Let  U  denote  the  set  of  allowed  control  values  in  R™.  Define  the  equilibrium  manifold 

M  =  {(x,u)  I  f(x,u)  =  0.  u  in  U,  and  f»(x,u)  has  rank  n). 

Projecting  M  onto  the  first  factor  we  obtain  M,  the  set  of  equilibrium  states.  Note  that  M  and  M  depend  on  the 
specified  control  limits. 

We  use  T  and  0  as  coordinates  on  M. 

Suppose  T  and  0  are  fixed.  Note  that  Q  =  0,  so  the  only  unspecified  terms  in  equation  (1.1)  are  the  aerodynamic 
forces  and  moment  However,  the  direct  aerodynamic  force  terms  in  the  U,W  degrees  of  freedom  must  exactly 
balatKe  the  gravity  and  thrust  forces  for  an  equilibrium  to  exist  ThCTefore,  the  aerodynamic  force  vector  in  the 
U.W  tangent  space  is  completely  determined,  even  though  we  do  not  know  the  point  (U,W). 

We  introduce  the  notion  of  residualized  aerodynamic  functions. 

Definition  2.1:  For  the  model  represented  by  equation  (1.1),  the  residualized  aerodynamic  functions  are 
C,(a),  Cj(a)  defined  by; 


.  . 

Cz(a) 

^  Ckw 

C*(a) 

p 

r 

(a) 


(2.2) 


The  meaning  of  the  residualized  aerodynamic  functions  is  simple:  at  any  equilibrium  state  the  aerodynamic  force 
vector  in  the  U.W  tangent  space  is: 


'f.(U.W) 

Fz(U.W) 


C,(a) 

Cz(a) 


(2.3) 


as  can  be  verified  by  solving  the  equation  Q  =  0  in  (1.1).  A  problem  may  arise  at  values  of  a  for  which  CM^(a) 
is  zero  or  (relatively)  very  small  in  magnitude.  We  deal  with  this  problem  by  splitting  it  into  two  subproblems: 
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1)  Asaime  unlimited  control  authority  and  work  with  the  full  system 

2)  Restrict  attention  to  a  subset  of  the  states  on  which  control  authority  is  adequate 

Denote  by  A  the  set  of  a  corresponding  to  equilibrium  states;  A  is  all  angles  from  0  to  27c  if  CM^(a)  never  van¬ 
ishes  and  5  is  unlimited.  In  practice,  elevator  deflection  angle  is  limited  and  its  control  effectiveness  degrades  at 
angles  of  attack  too  far  from  zero,  so  in  real  rqrplications  the  set  A  turns  out  to  be  a  proper  subset 

We  define  the  notion  of  univalent  residualized  aerodynamic  functions; 

Definition  2,2;  The  residualized  aerodynamic  functions  are  called  univalent  if  the  function  (Cx(a),Q(a))  defines 
a  one-to-one  mapping  from  A  onto  its  image  in  the  space  of  directions  in  the  U,W  tangent  space. 

Observation  2.1;  Suppose  the  residualized  aerodynamics  are  univalent  Then  for  each  fixed  value  of  T  and  6 
there  is  at  most  one  equilibrium  state  x(T,d)  in  M. 

The  validity  of  this  observation  follows  easily  firom  the  earlier  discussion.  Once  the  direction  of  the  aao- 
dynamic  force  at  the  equilibrium  is  known  then  the  value  of  a  is  uniquely  detomined.  Then  is  determined 
by  the  force  magnitude,  and  the  equilibrium  values  of  U  and  W  are  uniquely  determined. 

Section  3;  Dynamic  Inversion 

The  control  problem  we  consider  here  is  the  foUowing; 

Statement  of  the  problem;  Give  an  equilibrium  state  x,  determine  a  controller  u  =  K(x)  so  that  x  is  a  global 
attractor  for  the  system 

X  =  f(xaC{x))  (3.1) 

It  is  clear  that  any  global  attractor  for  (3.1)  must  be  an  equilibrium  state.  Using  dynamic  inversion  we  will 
address  this  problem  for  vehicle  models  having  univalent  residualized  aerodynamic  functions  (see  Section  2  for 
definitions). 

Dynamic  inversion  is  a  very  simple  technique  (sec  [E1],[E2],  and  [MEHH]).  We  assume  here  that  the  reader 
knows  the  basics. 

The  approach  we  take  is  to  invert  the  rotational  degrees  of  freedom  to  a  set  of  stable,  second  order  dynamics. 
We  leave  the  value  of  T  fixed  --  we  could  use  T  as  a  control  input  but  we  choose  not  to.  The  throttle  is  typi¬ 
cally  used  as  a  low-bandwidth  control  that  is  not  changed  during  dynamic  maneuvers.  In  the  rest  of  this  section 
we  show  how  K  is  constructed  and  discuss  stability  of  the  attitude  dynamics.  The  big  question  concerning  the 
stability  of  the  complementary  dynamics  is  addressed  in  Section  4. 

To  construct  the  controller  K,  first  select  a  desired  stable  set  of  second-order  linear  dynamics  for  6; 

Q  =  -  2  C®  Q  -  ©2  (0  _  0^  (32) 

Desirable  sets  of  d3mamics  of  this  type  for  aircraft  pitch  axis  control  are  described  in  the  MIL-F-878SC 
specifications  document  These  dynamics  ate  realized  in  the  closed-loop  system  if  5  takes  values  according  to 
the  following  equation; 
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(3J) 


*  I 

Qand  -  Ctm 
Cm^Co) 


where  Qaoa  the  value: 


=  -  2  Cw  Q  -  to^  (0  -  Oond) 


(3.4) 


and  6 _ _  is  the  6  value  of  the  desired  equilibrium  state  x.  These  equations  will  do  the  job  so  long  as  they  give 

allowed  values  of  5. 

The  surface  saturation  problem  at  high  a  appears  in  all  real  aircraft  That  is  the  motivation  for  using  thrust¬ 
vectoring  for  high-a  controL 

The  problem  of  attitude  control  by  dynamic  invasion  using  thrust  vectoring  was  solved  by  Mike  Elgersma  in 
[El].  For  our  current  problem,  we  assume  the  existence  of  an  invariant  set  in  the  closed-loop  state  space  inside 
which  the  the  pitch  dynamics  have  been  stabilized  by  dynamic  inversion  using  aerodynamic  surfaces  o  ly. 

Recoit  analysis  suggests  there  is  an  invariant  sa  closely  tied  to  practical  aircraft  flight.  This  set  is  a  bounded 
region  such  that: 

1)  a  lies  between  the  zero-lift  value  Oq  <  0  degrees  and  the  maximum  lift-to-drag  value  a„,„, 

2)  V  lies  between  a  minimum  cruising  speed  and  a  maximum  speed  V„... 

3)  the  flight-path  angle  y  s  0  -  a  lies  in  a  fixed  range  about  zero  (say  ±  10  degrees), 

4)  T  is  greater  than  or  equal  to  the  value  needed  to  fly  maximum  y  at  the  V^iJcu,..  condition 

5)  0  is  constrained  to  lie  in  an  interval  about  zero  consistent  with  conditions  1  and  3 

6)  Q  is  constrained  to  have  magnitude  smalla  than  a  suitable  e 

Section  4:  The  Stability  Result 

For  pitch-axis  models,  computer  simulatiois  show  that  the  pitch-attitude  control  strategy  presented 
has  excellent  stability  properties  ova  extreme  ranges  of  initial  and  transient  state  conditions.  So 
applied  the  technique  without  a  rigorous  proof  that  it  should  work  globally.  In  this  section  we  give 
biUty  result  that  begins  to  explain  the  results  observed  for  pitch-axis  applications. 

We  sketch  the  proofs  for  the  unlimited  control  authority  case. 

Lemma  4.1:  Suppose  the  aerodynamic  drag  is  always  positive.  Then  thoe  is  a  finite  disc  D  centered  at  the  ori¬ 
gin  in  the  U,W  plane  into  which  all  trajectories  eventually  enta  and  remain. 

Proof:  From  equation  (1.1)  we  can  compute  the  time  rate  of  change  of  -i-  W^: 

^  (U®  +  W^)  =  2  gVsin(a  -  0)  +  2  TVcos(a)  -  pV^S  CD(a,5).  (4.1) 

dt 

The  expression  Co  denotes  the  total  drag  coefficient  including  direct  surface  effects.  The  first  two  terms  on  the 
right  side  are  linear  in  V,  while  the  last  is  cuImc.  For  sufficiently  large  V  the  last  term  will  dominate. 

Lemma  A2:  La  F(U,W)  =  (F,(U,W),  F,(U,W))  denote  the  residualized  aerodynamic  force  vector  defined  in 
equation  2.3 

Define  the  residualized  system: 


in  Section  3 
far  we  have 
a  global  sta- 
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(4.2) 


V 

u 
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div(F)  = 


au 


(4.3) 


for  all  (U,W)  not  equal  to  (0.0).  Then  the  only  possible  closed  orbits  of  the  residualized  system  are  equilibria. 

Proof:  The  proof  is  by  contradiction.  Let  d>  denote  the  right  hand  side  of  equation  4.2.  Suppose  there  is  some 
closed  orbit  C  -  either  a  limit  cycle  or  a  homoclinic  orbit.  Let  A  denote  the  interior  of  C.  By  Green’s  theorem 
in  the  plane. 


J  c  «b.v>  =  J  A  div(<I>)  dArea 


(4.4) 


where  v  denotes  the  outward  pointing  unit  normal  vectt^  to  C.  The  integrand  on  the  left  vanishes  identically  by 
construction,  while  the  integrand  on  the  right  is  equal  to  div(F)  which  is  negative  by  assumption.  This  contradic¬ 
tion  proves  the  lemma. 

Remark  4.1:  This  Lemma  (and  prooQ  aie  given  as  exercise  1.3  in  chapter  2  of  [H].  We  approach  the  fxoblem 
this  way  because  some  examples  we  have  looked  at  (the  F-4,  F-14,  and  F-15  aircraft)  appear  to  satisfy  the  nega¬ 
tive  divergence  condition.  See  Figures  2  and  3  (copied  from  [MEHH]). 

We  have  already  observed  in  Section  2  (Ob%rvMion  2.1)  that  for  models  having  univalent  residualized  aero¬ 
dynamics  there  is  a  unique  equilibrium  state  when  T  and  6  are  fixed.  Under  the  additional  condiuon  that  the 
residualized  aerodynamic  vector  field  has  negative  divergence,  it  follows  by  the  Poincaire-Bendixson  theory  [H] 
that  when  T  and  6  are  held  fixed  and  the  aoodynamics  have  negative  divergence,  the  unique  equilibrium  state  is 
a  global  atuactor  for  the  U.W  dynamics  inside  the  invariant  set 

Section  5:  Linear  Models 

In  the  linear  case  the  issue  of  stability  is  related  to  right  half  plane  zeros  in  the  elevator  to  pitch  angle  transfer 
function.  Under  the  usual  assumptions  [McAG]  this  transfa  function  is  given  by 


A(s+^)(s-h^) 

(9^ _ Ti  Tz 

8  (s^+aj*s-hbph)(s?+a^+b^) 


(5.1) 


where  both  zeros  are  typically  in  the  left  half  plane  with  a  few  excq}tions:  there  are  a  few  conditions  whoe 

M 

is  negative  and  tterc  are  post  stall  conditions  where  —  is  negative. 

Ti 

Tte  application  of  dynamic  inversion  to  the  pitch  axis  in  the  linear  case  gives  a  closed  loop  characteristic  equa¬ 
tion  given  by 


♦d(s)=(s*+2CaB+<o*Ks+4-Ks+-:p) 

1]  I2 


(5.2) 
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so  right  half  plane  zeros  of  give  unstable  closed  loop  poles  but  yet 

o 


_0 _ ca^ 

0e  (s^+2Ctos+<o^ 


(5.3) 


where  to  and  ^  are  specified.  We  conclude: 

Observation  S.l:  the  dynamic  inversion  controllCT  discussed  in  Section  3  for  the  linear  model  is  stable  if  and 
only  if  T,  and  T2  are  positive. 

Observation  S.l  makes  it  easy  to  check  the  stability  of  the  closed-loop  system  by  looking  at  standard  tabular 
data  for  the  open-loop  zeros.  The  open-loop  parameters  fcv  the  transfer  functions  of  10  aircraft  at  94  total  flight 
conditions  is  tabulated  in  Heffley  and  Jewell’s  data  book  [HI].  Of  these  94  flight  conditions  there  were  four  with 

right  hall  piane  <.aos:  the  four  excqitional  values  were  in  die  range  -0.027<7^<-0.00049  rad/sec.  None  of  the 

M 

flight  conditions  were  post  stall  so  —  was  positive  for  all  94  conditions. 

‘2 

These  are  low  frequency  instabilities  that  pilots  or  autopilots  usually  stabilize  with  a  combination  of  additional 
outer  loop  feedbacks:  the  outer-loop  feedback  maps  a  combination  of  altitude,  altitude  rate,  and  velocity  to  a 
combination  of  elevator  and  throtde. 

In  [McAG]  there  are  ^proximate  fonnulas  for  that  show  that  it  is  the  variation  of  thrust  minus  drag  with 

M 

velocity  that  causes  the  right  half  plane  zeros  for  3  of  the  4  exceptional  flight  conditions.  Since  drag  increases 
with  velocity  (stabilizing)  we  conclude  the  right  half  plane  zeros  result  in  part  from  propulsion  chaim:teristics 
(coisequently,  the  stability  result  of  Section  4  should  not  be  expected  to  apply).  These  3  cases  were  F-4  at  Mach 
0.6  and  35K  ft  altitude,  and  the  XB-70  at  Mach  0.6  and  20K  ft  and  Mach  0.9  and  40K  ft  altitude. 

The  one  other  exceptional  case  is  intmsting.  That  case  is  die  X-IS  at  Mach  1.6  and  80K  ft  altitude.  The  right 
half  plane  zoo  is  not  associated  with  propulsion  because  the  engine  was  off.  We  do  not  know  why  the  stability 
result  of  Section  4  does  not  work  here  (other  than  to  say  that  some  hypothesis  was  violated).  This  flight  condi¬ 
tion  was  exceptional  in  that  the  angle  of  attack  above  14  degrees  wa:,  the  largest  of  all  those  listed  for  the  X-IS 
(10  total). 
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